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SUMMARY 

This  Memorandum  describes  a  computer  program  which  can  be  used  to  predict  the 
growth  rates  of  fatigue  cracks  emanating  from  notches  and  growing  through  residual  stress 
fields.  The  residual  stress  distributions,  alternating  loading  conditions  and  specimen 
geometry  must  be  specified  by  the  user.  The  program  uses  a  Green's  function  technique  to 
calculate  the  stress  intensity  factor  due  to  the  applied  and  residual  stresses  for  any  crack 
length.  The  calculated  stress  intensity  factor  can  be  corrected  to  account  for  the  exact  crack 
shape,  if  it  is  known.  The  crack  growth  rate  is  obtained  from  a  database  of  experimentally 
determined  crack  growth  data  as  a  function  of  stress  intensity  factor,  for  a  number  of 
different  materials. 
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INTRODUCTION 


The  most  common  sites  of  fatigue  damage  in  aircraft  structures  are  fastener  holes. 
The  use  of  life  enhancement  techniques  such  as  cold  expansion  and  interference  fit  fasteners 
are  now  commonplace  in  order  to  combat  this  problem.  Cold  expansion  increases  the 
fatigue  endurance  of  components  by  reducing  the  rate  at  which  fatigue  cracks  grow.  This  is 
achieved  by  inducing  a  compressive  residual  stress  field  around  holes  in  the  region  through 
which  the  fatigue  cracks  must  grow.  The  extent  of  the  compressive  residual  stress  zone  is 
typically  only  a  few  millimetres  and  outside  this  resides  a  balancing  tensile  residual  stress 
zone.  Most  of  the  fatigue  life  of  a  component  however  is  consumed  in  the  growth  of  small 
cracks  and  therefore  reducing  the  growth  rate  of  such  cracks  can  result  in  significant 
increases  in  overall  fatigue  endurance. 

The  design  of  new  aircraft  structures  or  repair  schemes  is  now  frequently  performed 
using  a  damage  tolerant  approach.  Initial  defects  of  specific  sizes  are  assumed  to  exist  at  all 
stress  concentrations  such  as  fastener  holes.  Fracture  mechanics  analyses  are  used  to 
calculate  the  fatigue  life  consumed  in  growing  a  crack  of  the  assumed  initial  defect  size  to  a 
size  which  can  be  reliably  detected  by  inspection.  Damage  tolerant  design  applied  to  cold 
expanded  holes  therefore  requires  that  the  growth  rates  of  cracks  subjected  to  both  residual 
and  applied  stresses  can  be  calculated. 

This  Memorandum  describes  a  method  of  calculating  the  growth  rates  of  cracks 
subjected  to  cyclic  loading  and  growing  from  open  holes  through  arbitrary  residual  stress 
fields.  It  also  describes  a  computer  program  which  can  be  used  to  perform  these 
calculations  from  information  input  by  the  operator.  A  Green's  function  technique  is  used 
to  calculate  stress  intensity  factors  at  crack  lengths  specified  by  the  user.  Residual  and 
applied  stresses  and  specimen  geometry  also  need  to  be  input  by  the  user.  The  total  stress 
field  is  calculated  by  summing  the  residual  and  applied  stresses  and  stress  intensity  factors 
are  calculated  by  integrating  the  product  of  the  stress  distribution  and  the  Green’s  function 
over  the  crack  length.  The  Green’s  function,  which  is  described  in  Ref  1 ,  is  appropriate  to  a 
through-the-thickness  crack  in  an  infinite  sheet;  corrections  need  to  be  made  for  the  more 
commonly  occurring  comer  crack.  Various  corrections  for  crack  and  specimen  geometry 
are  available  in  the  program  and  are  described  in  this  Memorandum.  Crack  growth  rates  are 
calculated  from  a  materials  database  for  BS  L65  and  BS  L71  aluminium  alloys  which  is 
included  in  the  program  and  can  be  added  to  by  simple  modifications.  The  method  of 
operation  of  the  computer  program  is  described  and  the  results  of  a  test  case  presented. 
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2  THEORETICAL  BACKGROUND 

2 . 1  The  prediction  method 

This  section  outlines  the  method  used  to  obtain  crack  growth  rates  from  a  knowledge 
of  the  residual  and  applied  stresses.  Residual  stress  distributions  induced  prior  to  fatigue 
loading  need  to  be  determined  and  stored  in  the  computer.  Section  2.2  describes  the  main 
features  of  residual  stress  distributions,  how  they  are  specified  and  how  the  program 
interpolates  values  at  any  location.  The  applied  stress  distribution  is  calculated  by  the 
program  as  described  in  section  2.3.  The  total  stress  distribution  (residual  plus  applied)  can 
thus  be  derived  and  stress  intensity  factors  calculated  with  the  use  of  an  appropriate  Green's 
function.  In  section  2.4  a  specific  Green’s  function  is  described  which  is  appropriate  to  a 
through  crack  growing  from  a  central  notch  in  an  infinite  component.  This  Green's  function 
is  used  in  conjunction  with  the  residual  and  applied  stresses  to  give  the  required  stress 
intensity  factor  (K)  solution.  The  numerical  integration  technique  used  to  give  the  K 
solution  is  described  in  section  2.5.  Since  the  Green's  function  is  appropriate  to  a  through 
crack  in  an  infinite  sheet,  the  K  solution  will  in  general  need  correcting  to  account  for  crack 
and  component  geometries.  The  corrections  available  in  the  program  are  described  in 
section  2.6.  Having  obtained  the  required  K  solution  it  only  remains  to  calculate  the 
appropriate  crack  growth  rates  at  specified  crack  lengths.  The  method  of  calculation  is 
described  in  section  2.7  along  with  a  description  of  the  database  used. 

2.2  Residual  stress  fields  and  their  interpolation 

The  residual  stress  fields  considered  in  this  Memorandum  are  appropriate  to  typical 
stress  fields  induced  around  open  holes  by  prestressing  or  cold  expansion.  The  fields 
considered  were  obtained  from  a  number  of  sources  such  as  finite  element  analyses, 
mathematical  models  describing  the  stress  fields,  or  simple  linear  approximations. 
Whichever  source  was  used,  the  distributions  have  certain  features  in  common  as  illustrated 
in  Fig  la.  These  are: 

(a)  a  compressive  residual  stress  at  the  hole  edge, 

(b)  the  residual  stress  decreases  in  magnitude  with  distance  away  from  the  hole 
becoming  tensile  and  reaching  a  maximum  a  short  distance  from  the  hole 
edge, 

(c)  the  tensile  residual  stress  then  decreases  with  distance  from  the  hole, 

(d)  the  magnitude  of  the  maximum  compressive  stress  is  significantly  greater 
than  the  magnitude  of  the  maximum  tensile  stress. 

Some  of  the  distributions  (see  Fig  lb)  also  display  a  region  referred  to  as  a 
compressive  yield  zone  near  the  hole  edge  where  the  compressive  stress  initially  increases  in 
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magnitude  with  distance  from  the  hole.  This  occurs  if  the  compressive  stress  at  the  hole 
exceeds  the  compressive  yield  point  of  the  material. 

Residual  stress  fields  are  usually  described  as  a  .set  of  values  at  fixed  coordinate 
points.  It  is  necessary  to  employ  an  interpolation  routine  to  calculate  values  between  those 
specified  and  a  Lagrangian  technique  suitable  for  use  on  a  computer  was  selected. 
Unfortunately  the  interpolation  was  not  sufficiently  accurate  at  "cenain”  points  in  the 
distribution,  for  example  at  sharp  changes  of  gradient  or  points  of  inflection,  so  that  the 
residual  stress  distribution  had  to  be  split  into  regions.  The  typical  residual  stress 
distribution,  shown  in  Fig  la,  is  split  into  three  regions  which  enables  accurate 
interpolation.  The  residual  stress  distribution  in  which  compressive  yielding  has  occurred, 
shown  in  Fig  I  b,  must  be  split  into  at  least  four  regions. 

2.3  Stress  intensity  factors 

A  stress  intensity  factor  is  a  single  parameter  which  describes  the  magnitude  of  the 
stress  field  in  the  vicinity  of  a  crack  tip.  The  assumption  is  usually  made  that  the  behaviour 
of  a  crack  is  controlled  by  the  stresses  around  its  tip,  hence  stress  intensity  factors  are  verj- 
useful  in  characterising  crack  behaviour  in  an  elastic  body  and  are  widely  used  for  predicting 
the  rate  of  growth  of  fatigue  cracks.  The  stressing  mode  appropriate  to  a  stress  intensity 
factor  is  usually  denoted  by  Kj^  ,  where  N  =  I ,  II  or  III  specifies  the  direction  of  relative 
movement  of  the  crack  faces  as  opening,  sliding  or  tearing  respectively.  The  cracks  with 
which  this  report  is  concerned  are  opening  mode  cracks  (N  =  I)  growing  along  a  line 
perpendicular  to  the  remotely  applied  stress. 

The  stress  intensity  factor  for  the  situation  of  a  through  crack  of  length  L  measured 
from  the  edge  of  a  circular  hole  of  radius  R  in  an  infinite  sheet  (see  Fig  2)  is  given  by: 

L 

Kj  =  l/VnLja(x)G(L/R,x/L).dx  ,  (1) 

0 

where  a(x)  is  the  stress  distribution  along  the  crack  site  in  the  uncracked  configuration, 

G(L/R,  x/L)  is  the  appropriate  Green's  function  for  this  configuration. 

The  stress  distribution  along  a  line  from  a  hole  in  an  infinite  sheet  perpendicular  to  a 
remotely  applied  stress  is  given  by: 


where  R  =  radius  of  the  hole 

X  =  distance  along  the  crack 
S  =  remotely  applied  stress 
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The  stress  intensity  factor  distribution  for  this  stress  field  is  known.  Values  of 
stress  intensity  factor  were  calculated  with  no  residual  stress  field  present  and  compared 
with  the  known  solution  to  test  the  accuracy  of  the  computer  program.  If  a  residual  stress 
field  is  present  in  the  component  then  the  stress  distribution  is  given  by  the  sum  of  the 
applied  and  residual  stress  distributions. 

2.4  The  Green's  function 

The  Green’s  function  used  was  developed  by  Shivakumar  and  Forman  ^  and  is  for  a 
through  crack  at  the  edge  of  a  circular  hole  in  an  infinite  sheet. 

The  Green's  function  G(L/R,  x/L)  can  be  expressed  as  G(A.,P) ,  that  is, 


’/2 

G(?i.P)  =  +  C(a,P) 


(1+F  ) 

2 - L 

(l-p¥ 


(— ) 

'^1  -B^ 


(3) 


where  X  =  L/R  ,  p  =  x/L  and  a  =  1/(1  +  X),  and 


and 


Fp  =  (1  -  p2)(.2945  -  .3912p2  +  .7635p4  -  .9942p6  +  .5094p8) 

/ 


C(a.  p)  =  Z 


m=l 


mil  an/2 


n=0 


where  are  constants  whose  vailues  are  given  in  Table  1. 


(4) 

(5) 


Table  1 

Coefficients  in  Green's  function 


m/n 

0 

1 

2 

3 

1 

0.8164 

-4.5911 

7.6059 

-3.8529 

2 

0.0492 

17.3181 

-36.8465 

20.6753 

3 

-0.4831 

-30.5563 

75.4833 

^.3540 

4 

-0.1746 

26.2877 

-73.1167 

44.2607 

5 

0.7952 

-8.4570 

26.8666 

-16.7296 

2.5  Numerical  integration 

It  can  be  seen  that  the  Green’s  function  is  algebraically  complex  and  for  this  reason, 
as  well  as  the  fact  that  the  residual  suc:»  distribution  is  unlikely  to  be  represented  by  a 
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simple  polynomial,  the  integration  has  to  be  carried  out  by  numerical  means.  The  meihcKi 
selected  is  a  Gaussian  quadrature  numerical  integration  method,  chosen  because  not  only  is 


it  very  accurate  but  also  because,  using  suitable  weighting  functions,  it  is  able  to  take 


account  of  the  singularity  which  occurs,  due  to  the  Vl  -  p  and  s/ 1  -  denominators,  in 


the  Green's  function  as  (i  1  . 


Gaussian  integration  is  carriea  out  as  follows:-  Given  a  weighting  function  VVT  Y) 
and  a  function  F(Y) ,  there  exists  a  Gaussian  formula  of  the  form 
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JW(Y)F(Y)dY=:  Ih..F(Y.)  , 


where  Y^  are  the  zeros  of  the  appropriate  nth  degree  orthogonal  polynomial,  and  h^  are  the 
corresponding  quadrature  coefficients.  Both  Y.  and  hj  are  tabulated  for  various  weighting 
functions^. 

For  the  purposes  of  this  work,  the  function  F(Y)  is  the  product  of  the  reduced 
Green’s  function  g(k,  P)  and  the  total  stress  distribution,  that  is  o(x)g(A,  p)  ,  and  the 
weighting  function  is  1/V(1  -  p),  where  g(X,  p)  is  defined  as 

g(X,  P)  =  GO.,  p)  (7) 

An  eighth  order  integration  (n  =  8)  was  chosen,  as  it  gave  an  acceptable  compromise 
between  running  time  and  accuracy.  Obviously,  the  greater  the  number  of  summations  that 
are  earned  out,  the  longer  will  be  the  processing  nme  but  also  the  greater  the  accuracy  of  the 
process.  However,  it  was  found  that  no  significant  increase  in  accuracy  was  attained  by 
carrying  out  a  higher  order  integration  than  eight,  so  the  increase  in  processing  time  required 
to  caiT)'  out  such  a  process  was  thought  to  be  unnecessary. 

Because  of  the  shape  of  some  of  the  residual  stress  fields  under  consideration,  it  was 
necessary  to  carry  out  the  integration  in  regions,  as  described  in  section  2.1.  This  is  due  to 
the  fact  that  the  ir.terpolation  method  was  unable  to  take  sufficient  account  of  sharp  changes 
in  gradient  of  the  stress  distribution;  the  distribution  was  therefore  interpolated  in  sections. 

To  this  end  the  following  scheme  was  worked  out:- 

The  integral  expression  for  Kj  in  equation  (1)  can  be  written 
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The  distribution  is  split  into  n  regions  and  the  boundaries  between  regions  are 
denoted  by  Lj,  L2,  with  =  L^/L  for  j  =  1.  2  ...{n  -  1);  =  0  and  p^  =  1 

The  width  of  a  region  is  given  by  =  P^  -  P(j_i)  •  Thus,  in  general. 


K 


I 


(9) 


where 


j  ^  I  g(L.p)g(X.p) 


vl  -  p 


(10! 


Substitution  of  y  =  (P  -  p^  jVA^  leads  to 


1  =  \ 


-  '  f^lL.yA^  +  P^_j)g(A„yA^  +  P^_j) 

1  -  yA .  -  p.  , 

^  I  J-i 


i 


dy  .  j  =  1 ,2,...(n  -  1 ) 


A. 

J 


and 


(11) 


I  I  o|L(.vA^-P„_|)lg(yA^->|i^  ^ 

"  ■  "0  '  ' 


The  numerical  integrations  for  Ij(j  =  l,...n-l)  were  done  using  Gaussian 
quadrature  with  unit  weighting  function  as  the  integrand  is  not  singular.  Integration  of 
was  done  using  Gaussian  quadrature  with  a  weighting  function  of  (1  -  y)“  •  to  allow  for 
the  singularity  at  y  =  1  . 

2.6  Corner  crack  corrections 

As  has  already  been  observed,  the  stress  intensity  factors  obtained  by  the  Green's 
function  approach,  are  for  through  cracks  in  an  infinite  sheet.  In  order  to  obtain  stress 
intensity  factors  for  the  more  commonly  occuring  comer  cracks  in  a  sheet  of  finite 
dimensions,  correction  factors  need  to  be  applied.  The  correction  factors  relate  to  the  shape 
of  the  crack  and  the  width  and  thickness  of  the  sample  in  which  the  crack  is  growing.  Both 
correction  factors  (crack  shape  and  finite  width)  are  based  on  the  work  of  Newman  and 
Raju^  who  developed  empirical  equations  for  stress  intensity  factors  based  on  finite  element 
analyses.  The  equations  developed  are  relevant  to  quarter-elliptical  comer  cracks  in  finite 
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width  compone-.  s.  They  account  for  both  component  width  and  thickness,  and  the  shape 
of  the  crack  ’  .  terms  of  its  aspect  ratio  (bore  length/surface  length).  Since  the  Newman  and 
Raju  worK  is  for  the  case  of  no  residual  stress,  then  the  correction  factors  required  cannot  be 
directly  determined.  Two  methods  are  used  to  determine  the  comer  crack  shape  corrections 
and  one  method  to  determine  the  finite  width  correction. 

Firstly  examining  the  factors  in  equations  (50-63)  of  Ref  3.  it  is  apparent  that  the  g-, 
factor  in  Ref  3  accounts  for  the  presence  of  the  hole  in  calculating  the  stress  intensity  factor. 
It  is  very  similar  in  magnitude  to  the  factor  derived  using  the  Shivakumar  and  Forman 
method  assuming  no  residual  stress.  The  first  correction  factor  offered  in  the  present 
program  is  the  solution  calculated  by  Newman  and  Raju  with  the  term  set  to  1 ,  This 
correction  factor  is  multiplied  by  the  Shivakumar  and  Forman  solution  for  a  through  crack  in 
a  residual  stress  field.  This  gives  the  solution  for  an  elliptical  crack  in  an  arbitrary'  residual 
stress  field  in  an  infinite  width  component.  This  will  subsequently  be  referred  to  as  the  G2 
correction.  The  remaining  terms  in  the  equations  refer  to  the  geometnes  of  the  test  piece  and 
crack.  The  crack  aspect  ratio  and  the  specimen  thickness  have  to  be  input  by  the  user.  This 
method  assumes  that  the  g^  term  entirely  accounts  for  the  presence  of  the  hole  and  that 
there  is  no  'fine  tuning'  included  in  the  other  terms  in  order  to  achieve  a  gcxxi  fit  to  the  finite 
element  solutions. 

The  second  comer  crack  correction  factor  is  obtained  using  the  complete  Newman 
and  Raju  solution,  ie  including  the  g2  term.  The  solution  again  requires  the  input  of  both 
crack  aspect  ratio  and  component  thickness,  by  the  user.  The  complete  Newman  and  Raju 
solution  is  divided  by  the  Shivakumar  and  Forman  solution  assuming  no  residual  stress 
field  is  present.  This  results  in  a  correction  factor  which  when  multiplied  by  the 
Shivakumar  and  Forman  solution  with  a  residual  stress  field  present,  gives  a  solution  for  an 
elliptical  crack  in  a  residual  stress  field  in  a  component  of  infinite  width. 

There  is  also  an  option  in  the  program  to  correct  for  the  finite  width  of  the 
component.  The  correction  factor  available  is  that  due  to  Newman  and  Raju  and  is  given  by 
equation  (47)  of  Ref  3. 

2.7  Crack  rate  predictions 

Crack  rates  are  calculated  by  comparing  the  derived  stress  intensity  factors  with 
known  crack  rate  data.  The  program  calculates  the  maximum  and  minimum  stress  intensity 
factor  values  and  hence  the  values  of  AK  and  stress  intensity  factor  ratio  R  .  A  subroutine 
interpolates  or  extrapolates  as  appropriate  from  a  database  of  crack  rates  versus  stress 
intensity  factors  at  a  range  of  R  values  fora  specific  material  to  calculate  crack  growth  rate. 
The  routine  used  in  the  program  is  that  developed  by  Edwards'^.  The  database  used  is  also 
the  one  described  in  Ret  4  and  k  based  on  data  obtained  using  BS  L71  and  BS  L65 
aluminium  alloys.  The  data  is  con.ained  in  Lines  10150  to  10290  in  the  program.  Other 
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data  may  be  substituted  in  the  correct  form  (see  Ref  4).  or  the  program  mcxisfied  lo  access  an 
alternative  database  of  CTack  growth  rates. 

3  THE  PROGRAM 

3.1  Program  outline 

The  program  was  written  in  Hewlett  Packard  BASIC  (BASIC  2.0)  and  i.s  listed  in 
the  -Appendix,  An  outline  of  the  scope  of  the  program  is  described  in  this  section.  Details 
of  how  to  run  the  program  are  given  in  section  3.2  which  can  be  used  as  a  users  guide.  The 
file  format  required  for  the  residual  stress  data  is  described  in  section  3.3.  This  foratat  must 
be  used  if  the  user  creates  data  files  which  can  be  accessed  by  the  program. 

The  program  will  calculate  stress  intensity  factors  at  crack  lengths  specified  by  the 
operator,  taking  into  account  a  residual  stress  field  if  required.  When  running  the  program, 
the  operator  will  need  to  input  the  following  information: 

(a)  Residual  stress  field;  the  residual  stress  field  can  be  defined  either  by 
inputing  coordinates  via  the  keyboard  or  by  loading  data  from  a  pre-srored  file,  the 
format  of  these  files  is  given  in  section  3.3.  It  is  also  possible  to  run  the  program 
with  no  residual  stress  field. 

(b)  Crack  lengths:  stress  intensity  factors  can  be  calcu.ated  either  over  a  range 
of  crack  lengths  at  regular  increments  or  for  a  number  of  crack  lengths  specified  by 
the  operator. 

(c)  Dimensions  of  the  component  and  crack  aspect  ratio:  stress  intensity  factors 
can  be  corrected  for  the  component  dimensions  (thickness  and  width)  and  the  crack 
aspect  ratio. 

(d)  Fracture  toughness  Kj^;  the  fracture  toughness  i.s  required  to  detennine 
when  failure  of  the  test  piece  will  occur. 

Stress  intensity  factor  values  will  be  calculated  at  the  specified  crack  lengths  and 
printed  out.  The  program  will  then  calculate  the  crack  growth  rates  corresponding  to  the.se 
stress  intensity  factors. 

The  operator  is  given  several  final  options  such  as  re-running  the  program  with  'he 
same  residual  stress  field,  storing  the  residual  stress  field,  stress  intensitj-  factors  and  rate 
values  for  future  reference,  plotting  current  data  and  plotting  previously  stored  distributions 
against  one  another  for  comparison. 

3.2  Running  the  program 

A  description  is  given  below  of  how  to  run  the  program,  including  all  of  tne  options. 
The  prompts  given  by  the  program  are  written  in  block  capitals  and  responses  in 
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parentheses.  The  running  procedure  is  split  into  six  sections  rhcse  sections  arc 

shown  in  the  flow  diagram  in  Fig  3. 

(a)  Inputing  residual  stres,s  fields  (RSF) 

(RUN) 

MULTI -PLOT  (Y/N)  -  the  operator  can  elect  to  run  the  routine  :o  plot  previously  stored 
distributions.  The  Y  response  is  described  at  the  end  of  this  section.  The  .N  response 
gives  the  following  prompt. 

INITIAL  OPTIONS  (Optl) 

0.  NO  RESIDUAL  STRESS  FIELD 

1.  INTERPOLATED  RESIDUAL  STRESS  FIELD 

2.  RESIDUAL  STRESS  FIELD  LOADED  FROM  FILE  (0.1  or  2) 

0.  to  run  program  without  any  residual  stress  field,  ie  with  applied  stress  only 

1.  to  input  a  residual  stress  field  from  the  keyboard, 

2.  to  input  a  residual  stress  field  from  a  previously  stored  file. 

Whichever  option  is  taken,  the  prc'gram  next  asks  for  the  applied  stress  range; 

MAX  REMOTE  STRESS  (MPai 

MIN  REMOTE  STRESS  (MPa) 

If  Optl  was  0,  then  the  prompt; 

RADIUS  OF  THE  HOLE  (mm)  is  given  and  the  operator  should  input  the  radius  of  the 
hole.  As  no  RSF  is  to  be  input,  the  program  goes  on  to  the  next  section  (*). 

If  Optl  was  1 ,  then  the  prompts  to  input  the  residual  stress  field  (RSF)  are  as  follows; 

RADIUS?  (mm)  -  operator  inputs  the  radius  of  the  hole, 

ORIGIN  OF  DISPLACEMENT  COORDS.  (OOR  1) 

0  -  FOR  CENTRE  OF  HOLE. 

1  -  FOR  EDGE 

VALUE  OF  X?  (mm) 

VALUE  OF  STRESS?  (MPa) . 


-  operator  inputs  the  origin  used  when  descnbing  the 

distances  in  the  residual  stress  distribution.  Jf 

-  operator  inputs  a  pair  of  coordinates  of  the  RSF;  the 
position  and  magnitude  of  the  residual  stres  . 

-  the  program  will  continue  to  give  these  prompts  after 
each  pair  of  coordinates  have  been  input  until  a  single 
carriage  return  is  input  lor  X. 
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The  coordinates  are  printed  on  the  screen  as  they  are  input  and  each  pair  given  a 
number.  When  the  input  ioop  is  terminated,  the  prompt; 

CHANGE  ANY?  (Y/N)  is  given.  The  operator  has  the  chance  to  change  any  pair  of 
coordinates  which  may  have  been  erroneously  entered.  If  there  are  entries  to  be  changed 
(Y),  the  prompt; 

ENTRY  TO  BE  CHANGED?  is  given.  The  operator  should  enter  the  number  of  the 
coordinate  pair  which  are  to  be  altered.  The  operator  will  then  be  asked  to  repeat  the  entrs 
number  and  should  enter  the  revised  coordinates.  The  operator  will  again  be  asked  if  their 
are  any  inputs  to  change  and  this  will  be  repeated  until  a  negative  response  (N)  is  given. 

STORE  DISTRIBUTION?  (Y/N).  The  operator  can  choose  to  store  the  RSF  in  a  tile  for 
future  use.  If  so; 

CREATE  NEW  FILE?  (Y/N).  The  RSF  can  be  stored  either  in  a  new  file  or  in  a  file  that 
already  exists. 

FILE  NAME?  The  operator  should  input  the  name  of  the  file  to  be  used/created. 

If  Optl  had  been  2  then  the  operator  would  be  asked; 

NAME  OF  FILE?  The  operator  inputs  the  name  of  the  file  containing  the  RSF  to  be  used. 

At  this  point,  the  Optl  =  1  and  Optl  =  2  sequences  become  coincident. 

(b)  Plotting  and  checking  of  residual  stress  fields 

The  coordinates  of  the  distribution  are  printed  on  the  screen.  Whenever  a  plot  is  to 
be  made,  from  any  part  of  the  program,  the  operator  will  be  given  the  range  of  coordinates 
on  both  axes  and  asked  to  define  the  size  of  the  intervals  within  these  ranges  at  which  to 
print  the  grid  dimensions. 

X-AXIS  LABEL  SPACING?  -  Operator  inputs  the  interval  at  which  grid  dimensions  on  X 
and  Y  axes  are  to  be  printed. 

Y-AXIS  LABEL  SPACING?  -  The  distribution  is  plotted,  after  which  the  operator  is  given 
the  chance  to  obtain  a  hard  copy. 

NO.  OF  REGIONS  -  The  operator  must  decide  on  the  number  of  regions  into  which  the 
distribution  is  to  be  split. 

INPUT  NUMBERS  OF  POSITIONS  OF  BOUNDARIES. 

POSITION  OF  END  OF  REGION  N  (where  N  =  No.  of  regions;  N  >  1 ) 

The  coordinates  of  the  RSF  values  are  printed  out  and  numbered  so  that  the  operator 
can  use  the  numbers  to  indicate  where  the  boundaries  between  the  regions  are  to  be  insened. 


The  distribution  is  now  re-plotted  this  time  using  data  interpolated  by  the  program. 
The  operator  is  thus  able  to  see  if  the  interpolated  data  accurately  represents  the  original 
input  data.  The  operator  is  asked; 

OK?  (Y/N)  and  should  indicate  whether  or  not  the  data  are  satisfactory. 

If  not,  the  program  allows  the  operator  to  re-position  the  boundaries  in  order  to 
produce  a  better  result.  The  program  therefore  returns  to  the  NO.  OF  REGIONS  prompt. 

(c)  Input  of  crack  lengths 

(*)  Now  that  the  RSF  has  been  input  (or  the  no  RSF  option  taken)  the  program  goes  on  to 
request  the  crack  lengths  at  which  stress  intensity  factors  are  to  be  calculated. 

1  SIFs  CALCULATED  OVER  A  RANGE  OF  CRACK  LENGTHS 

2  SIF  CALCULATED  FOR  SPECIHC  CRACK  LENGTHS 

If  option  1  (OPT2=l)  is  taken,  the  operator  will  be  asked  to  input  the  range  over 
which  the  SEFs  are  to  be  calculated,  thus; 

START  OF  RANGE,  (mm) 

END  OF  RANGE,  (mm) 

INCREMENTAL  LENGTH,  (mm) 

If  option  2.  is  taken  then  the  prompt; 

CRACK  LENGTH?  is  given  and  the  operator  can  input  one  or  a  list  of  crack  lengths.  The 
operator  must  input  a  crack  length  followed  by  a  carriage  return.  When  all  crack  lengths 
have  been  input,  the  response  to  the  CRACK  LENGTH  prompt  should  be  carriage  return. 

(d)  Correction  factors 

The  operator  is  now  asked  to  decide  on  the  type  of  correction  routine  to  be  used. 
The  options  are; 

0.  NO  CORRECTION 

1 .  NEWMAN  &  RAJU  WITHOUT  G2 

2.  N&R/SHRIV  AS  CORR 

For  an  explanation  of  the  two  corrections,  see  section  2.6. 

Unless  option  0.  is  taken,  the  dimensions  of  the  test  piece  will  be  required; 

THICKNESS?  (mm)  -  Operator  inputs  specimen  thickness 

A/C?  -  Operator  inputs  aspect  ratio  of  the  crack,  where  A  =  length  down  the  bore, 

C  =  length  along  the  surface. 


DO  YOU  WANT  A  WIDTH  CORRECTION?  (Y/N)  -  operator  is  given  the  option  of 

leaving  out  the  width  correction. 

WIDTH?  (mm)  -  if  required. 

(e)  Stress  intensity  factor  and  crack  growth  rate  calculations 

The  program  now  calculates  and  prints  out  the  maximum  and  minimum  stress 
intensity  factors  and  also  AK(=  K^ax  -  ^min  )  and  R(=  Kmm  /K^ax  )  After  each 
printout,  the  option  is  given  of  obtaining  a  hard  copy.  The  operator  has  the  option  of 
plotting  Kmax  vs.  L  and/or  AK  and  R  vs.  L  via  the  prompts. 

PLOT  KMAX?  (Y/N)  and 

PLOT  KMAX  -  KMIN  AND  R?  (Y/N) 

The  operator  is  then  asked  whether  rate  values  are  to  be  calculated; 

CALCULATE  RATE  VALUES?  (Y/N).  If  the  operator  responds  in  the  affirmative,  he  is 
given  the  prompt; 

KC  VALUE?  (MNm~^/2)  and  should  input  the  Kj^  value  for  the  material  under 
consideration.  The  program  then  calculates  and  prints  out  the  predicted  rate  values 
corresponding  to  the  calculated  stress  intensity  factors.  The  option  is  given  for  a  hard  copy 
of  the  rate  values  and  for  a  plot  of  rate  vs.  crack  length. 

The  operator  is  now  asked  whether  or  not  the  distributions  (RSF,  SIF's  and  rates) 
are  to  be  stored  for  future  reference  and/or  comparison  with  other  distributions.  If  the 
response  is  yes,  then  the  operator  can  either  store  them  in  an  existing  file  or  create  a  new 
one,  giving  the  name  of  the  file  to  be  used/created. 

STORE  DISTRIBUTIONS?  (Y/N) 

CREATE  NEW  FILE?  (Y/N) 

FILENAME? 

(f)  End  options 

ANOTHER  RUN  SAME  DIST?  (Y/N).  This  allows  the  program  to  be  re-run  with  the  same 
residual  stress  distribution,  for  example  if  different  crack  lengths  are  to  be  used.  If  the 
response  is  yes  (Y),  the  program  returns  to  (*). 

If  not  (N),  the  operator  is  now  given  the  prompt; 

PLOT  GROUPS  OF  DISTS?  (Y/N).  This  enables  the  operator  to  plot  one  or  more 
distributions  of  residual  stress  fields,  stress  intensity  factors  or  crack  growth  rates  from 
previously  stored  files. 

NAME  OF  FILE?  The  operator  should  enter  the  name  of  one  of  the  files  to  be  used. 
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ANOTHER  FILE?  (Y/N).  If  another  file  is  to  be  read  the  operator  enters  Y  and  the  program 
asks  for  a  file  name  again.  This  process  is  repieated  until  no  more  files  are  to  be  read  and  the 
operator  answers  N. 

These  options  are  then  given; 

o-Exrr 

1  -  RESIDUAL  STRESS  FIELD 

2  -  STRESS  INTENSITY  FACTOR  DIST. 

3  -  CRACK  GROWTH  RATE  DIST. 

The  operator  can  choose  which  of  the  three  distributions  to  plot  and  they  will  then  be 
plotted  on  the  same  axes  for  all  the  files  selected. 

Option  0.  ends  the  program. 

3.3  File  format 

Residual  stress  fields  may  be  input  via  the  keyboard,  as  described  in  section  3.2.  or 
may  be  accessed  by  the  program  from  a  data  file.  It  is  important  when  creating  a  data  file  to 
ensure  that  it  has  the  correct  format.  Details  of  the  file  format  required  are  given  below. 

The  general  form  of  the  data  file  is:- 

Lz  (1,0) 

Lz  (2,0) 

Lz  (3,0) 

Lz  (0,1),  Lz  (0,2) 

Lz(l,l),  Lz(l,2) 

Lz  (2,1),  Lz  (2,2) 

Lz  (n,l),  Lz  (n,2) 

Where  Lz  ( 1 ,0)  describes  the  coordinate  system  and  is  equal  to 

0  if  the  X  coordinates  are  measured  from  the  centre  of  the  hole 
1  if  the  X  coordinates  are  measured  from  the  edge  of  the  hole 
Lz  (2,0)  =  No  of  entries  in  the  file  array 
Lz  (3,0)  =  radius  of  the  hole 

Lz  (0,1)  to  Lz  (n,l)  are  the  distances  from  the  defined  origin 

Lz  (0,2)  to  Lz  (n,2)  are  the  residual  stresses  at  the  corresponding  positions 

The  maximum  dimension  of  the  array  is  (50,2),  if  more  than  50  coordinates  are 
used,  then  the  d  mension  statement  (line  7920)  must  be  altered. 


i 
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4  EXAMPLE  -  A  PRESTRESSED  HOLE 

An  example  of  the  use  of  this  program  is  presented  in  this  section  to  illustrate  its 
capabilities.  The  example  chosen  is  that  of  a  BS  2L65  aluminium  alloy  test  specimen 
containing  a  central  hole  (see  Fig  4).  Residual  stresses  may  be  induced  in  the  component  by 
applying  an  axial  load  large  enough  to  cause  local  yielding  around  the  central  hole.  A 
number  of  specimens  were  subjected  to  such  prestress  treatments  and  then  to  cyclic  fatigue 
loading  of  110  ±  96.5  MPa  ,  when  crack  growth  measurements  were  made  5.  The  crack 
growth  data  measured  in  the  experimental  investigation  will  be  compared  with  the  crack 
growth  data  predicted  using  this  computer  program.  The  specimen  geometry  and  loading 
conditions  which  need  to  be  input  to  the  program  are  therefore  defined  as  those  used  in  the 
experimental  investigation.  The  residual  stress  fields  resulting  from  various  prescress  levels 
were  determined  using  a  finite  element  analysis.  Other  inputs  required  by  the  program  will 
be  described  later  in  this  section. 

Three  values  of  prestress  were  used  in  the  experimental  investigation  resulting  in 
three  different  residual  stress  distributions  across  the  minimum  section  of  the  components. 
The  values  of  prestress  used  were  a  =  60%  of  the  yield  stress  on  the  net  section,  b  =  70% 
and  c  =  80%  ,  resulting  in  local  yielding  in  all  cases  at  the  edge  of  the  hole.  Residual 
stresses  determined  using  a  finite  element  program  are  plotted  in  Figs  5,  6  and  7  for  the 
three  prestresscs  ((a),  (b)  and  (c))  used  in  the  experimental  investigation.  Each  distribution 
was  split  into  three  or  more  regions  and  the  polynomial  expressions  derived  by  the  program 
are  also  shown  in  these  figures  as  solid  curves. 

Stress  intensity  factors  were  evaluated  at  crack  lengths  covering  two  ranges;  0. 1  to 
0.5  mm  in  steps  of  0. 1  mm  and  1  to  6  mm  in  steps  of  0.5  mm.  The  ranges  covered  crack 
lengths  at  which  experimental  crack  growth  data  had  been  measured.  Corrections  to  the 
stress  intensity  factor  were  made  using  both  of  the  comer  crack  routines  in  conjunction  with 
the  finite  width  correction  (described  in  section  2.6).  The  aspect  ratio  selected  was 
A/C  =  1 .4  which  was  thought  to  most  accurately  represent  what  occurred  in  practice.  Crack 
growth  rates  were  then  calculated  for  each  of  the  prestress  levels,  crack  length  ranges  and 
correction  factor  routines  described  above.  The  predicted  rates  are  shown  for  the  three 
prestress  levels  in  Figs  8,  9  and  10.  The  experimental  data  are  also  shown  as  solid  bars 
representing  the  range  of  measured  crack  growth  rates.  The  dots  on  each  of  the  solid  bars 
represent  the  log  mean  value  of  the  experimental  data. 

5  DISCUSSION 

The  residual  stress  distributions  input  to  the  program  consisted  of  eleven  coordinate 
values  and  are  shown  in  Figs  5  to  7.  The  polynomials  fitted  in  segments  by  the  computer 
program  were  a  reasonable  interpretation  of  the  input  data.  However  due  to  the  small 
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number  of  coordinates  specified,  abrupt  changes  in  slope  of  the  fitted  curves  occur  near  to 
the  region  boundaries  where  the  data  is  sparse.  It  is  expected  that  in  most  applications  a 
larger  number  of  residual  stress  coordinates  will  be  input  to  the  program  enabling  better  data 
fits.  Predicted  crack  growth  rates  using  the  two  different  correction  routines  (Figs  8  to  10) 
are  in  reasonable  agreement  with  each  other  and  with  the  experimentally  measured  data. 
This  gives  some  confidence  in  the  methods  used  but  does  not  give  an  indication  of  which 
correction  method  should  be  recommended. 

There  are  two  main  criticisms  of  the  approach  used  in  this  paper.  Firstly  it  is 
assumed  that  the  residual  stress  field  present  in  an  uncracked  body  remains  unaltered  by  the 
presence  and  growth  of  fatigue  cracks.  This  is  clearly  not  possible  and  the  residual  stress 
field  must  redistribute  as  cracking  occurs.  This  aspect  is  currently  being  investigated. 
Some  modification  to  the  models  used  in  the  program  may  be  possible,  dependent  on  the 
complexity  of  the  redistribution  process.  The  second  main  criticism  is  that  closure  of  the 
crack  flanks  is  not  taken  into  account.  Calculation  of  stress  intensity  factors  is  based  on  the 
range  of  the  total  local  stresses  and  takes  no  account  of  the  fact  that  the  crack  may  be  closed 
due  to  compressive  residual  stresses  acting  behind  the  crack  tip.  The  stress  intensity  factor 
range  experienced  by  the  crack  tip  will  be  different  from  that  predicted  using  the  methods 
described. 

Two  further  criticisms  of  the  method  are  also  being  considered.  The  first  is  that  the 
majority  of  the  fatigue  life  of  a  component  containing  beneficial  compressive  residual 
stresses  in  the  area  of  crack  initiation,  will  be  dominated  by  short  crack  growth.  It  has  been 
shown^  that  shon  cracks  propagate  faster  than  long  cracks  when  subjected  to  the  same 
stress  intensity  factor  range  calculated  by  conventional  linear  elastic  fracture  mechanics,  the 
effect  being  most  marked  at  negative  stress  ratios.  The  crack  rate  database  contained  in  the 
program  was  derived  from  long  crack  data  but  is  used  to  predict  the  growth  rates  of  short 
cracks.  From  the  foregoing  it  is  clear  that  predicted  growth  rates  at  short  crack  lengths  may 
be  underestimated  particularly  since  the  presence  of  compressive  residual  stresses  leads  to 
highly  negative  stress  ratios.  Shon  crack  growth  rate  databases  are  being  generated  for 
future  use  in  the  program.  The  second  criticism  is  that  the  Green's  function  solution  used  in 
the  program  is  not  applicable  to  most  engineering  situations  as  cracks  are  not  generally 
through  the  thickness  and  the  components  are  obviously  not  infinitely  wide.  The  corrections 
needed  to  modify  the  calculated  stress  intensity  factors  could  lead  to  considerable  errors  in 
the  predictions  if  they  are  not  entirely  appropriate.  Other  Green's  function  solutions  are 
currently  being  sought  and  determined. 


J 


18 


6  CONCLUSIONS 

(1)  A  method  has  been  devised,  and  a  computer  program  has  been  developed,  for 
predicting  the  growth  rates  of  fatigue  cracks  emanating  from  open  holes  and  propagating 
though  arbitrary  residual  stress  fields. 

(2)  The  program  has  been  used  to  predict  the  growth  rates  of  cracks  from  holes 
subjected  to  prestressing.  Reasonable  agreement  was  found  between  predicted  and 
experimental  results. 

(3)  Shortcomings  of  the  prediction  method  have  been  identified  and  suitable 
modifications  proposed. 
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Appendix 

PROGRAM  LISTING 


100  COM  /Opt/  Opt  I ,0pt2 

110  COM  /ftps/  Asmax  .flsnin 

120  COM  /Fie/  File(2  .50.2  ) 

130  COM  /Title/  Th  .B .Aoc  .Nanelt 30] 

140  DIM  6auss< 50  ,7 ) .Rpt( 50 .2  ) 

150  Th=0 
160  B=0 

170  Aoc=0 
180  K2=0 

1 30  Narne$= 

200  PRINTER  IS  1 
210  OUTPUT  2i"K"  i 
220  GClEAR 

230  CALL  Arrays 
240  PRINT 

250  PRINT 

251  01="" 

260  INPUT  "MULTI-PLOTP" ,0$ 

320  IF  0$="Y"  THEN 
330  CALL  Mplot 

340  GOTO  3020 

350  END  IF 

360  IF  0$="N"  THEN  GOTO  390 
370  GOTO  251 

33»/  PRINT 

400  ' 

410  I  INITIAL  OPTIONS 

420  PRINT  "OPTIONS;" 

430  PRINT 

440  PRINT  -0.  NO  RESIDUAL  STRESS  FIELD." 

450  PRINT 

450  PRINT  ’1.  INTERPOLATED  RESIDUAL  STRESS  FIELD" 

470  PRINT 

480  PRINT  "2.  RESIDUAL  STRESS  FIELD  LOADED  FROM  FILE- 

490  PRINT 

491  0S="" 

500  INPUT  0$ 

550  ON  ERROR  GOTO  491 

570  Optl=UAL(OS) 

580  OFF  ERROR 

590  IF  Opt1<'>0  AND  OptlOl  AND  0pt1''>2  THEN 
600  GOTO  491 

510  END  IF 

620  INPUT  "MAX.  REMOTE  STRESS'?  <MPa)".Max$ 

530  IF  LEN(Max$i=0  THEN  GOTO  520 

640  ON  ERROR  GOTO  620 

550  Asnax=UAL( Max$  ) 

560  OFF  ERROR 

570  INPUT  "MIN.  REMOTE  STRESS'?  <MPa)".MinS 
580  IF  L£N(Min$)=0  THEN  GOTO  670 

S90  ON  ERROR  GOTO  670 

700  Asnin=UAL( Min$  ) 

710  OFF  ERROR 

720  IF  Opt  1-0  THEN  SOTO  750 

730  CALL  Input s ( Radius  ) 

731  GCLEAR 

740  SOTO  810 

750  Nane$-"NO  RSF" 

750  INPUT  "RADIUS  OF  THE  HOLE'?  (nn)-,RadS 
770  IF  LEN(RadS)-0  THEN  GOTO  760 

780  ON  ERROR  GOTO  750 

790  Radius-UAL(Rad$ ) 

800  OFF  ERROR 

810  GCLEAR 
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820  OUTPUT  2; "K"  ; 

830  I 

840  'OPTIONS  ON  THE  TYPE  OF  RUN 
850  PRINT  "OPTIONS; • 

860  PRINT 

870  PRINT  "1.  SIP'S  CALCULATED  OUER  A  RANGE  OF  CRACK  LENGTHS. 
8S0  PRINT 

890  PRINT  "2.  SIF  CALCULATED  FOR  SPECIFIC  CRACK  LENGTHS." 

891  0S="" 

900  INPUT  0$ 

9S0  0pt2=VAL(0S) 

960  OFF  ERROR 

970  IF  OptCOI  AND  Opt2:.'2  THEN 

980  GOTO  891 

990  END  IF 

1010  I 

1020  I 

1030  CALL  Si f H Radius  .Gauss (•)  .Lp  ' 

1040  > 

1050  Th$=UAL$<Th> 

1080  B$=UAL$<B) 

1070  Aoc$=UAL$( Aoc > 

1080  1 

1090  IF  Th'0  THEN  Th$= " NC'l£ •' 

1100  IF  B=0  THEN  B$=“N0NE'' 

1110  IF  Aoc=0  Then  Aocs=“none" 

1120  I 

1130  ' 

1131  01=  ■■ 

1140  input  "DO  YOU  UANT  A  COPY’"  ,QJ 
1200  IF  0$='V-  THEN  GOTO  1240 
1210  IF  0S='N"  THEN  GOTO  1490 
1220  GOTO  M3I 
1230  I 

1240  PRINTER  IS  701 


1250 
1280 
1270 
1280 
1290 
1  300 
1310 
1320 
1330 
1  540 
1350 
1  360 
1370 
1  380 
1390 
1  400 
1410 
1420 


I 

I 

PRINT  NaiseS 

PRINT  "STRESS  INTENSITY  FACTOR  DISTRIBUTION." 

PRINT  DATES! TIMEDATE  ) 

PRINT  "RADIUS*  " t Radius :"nn," 

PRINT  "APPLIED  STRESS  RANGE*  " ; Asnax i jAsnin ; "MPa . " 

PRINT  "SPECIMEN  THICKNESS*  "iThS;"  hh." 

PRINT  "SPECIMEN  UIDTH=  "iBSi’  nn." 

PRINT  "A/C  RATIO*  " lAocS 

PRINT 

PRINT 

I 

PRINT  "CRACK  LENGTH"  , "KMAX . ( MNn" -3/2  )" ."KMIN"  ."KMAX-KMIN"  . "KM IN/KMAX " 
PRINT  "  ( nn )" 

PRINT 

FOR  Pr=0  TO  Lp 

PRINT  USING  "3X .DD.DO  .  14X  , SO. ODE  .7X  ,SD.2DE ,4X  .SD.2DE  .2X  . SD. ODD " i Gauss! Pr 


)  .GausstPr  ,2  )  ,6auss<Pr  ,3  )  .Gauss  1  Pr  ,4  )  .Gauss  (  Pr  .5  ) 


1430 

NEXT  Pr 

1440 

FOR  Iou=l 

TO  4 

1450 

PRINT 

1480 

NEXT  lou 

1470 

PRIMTEP  IS 

I 

1480 

1 

1490 

1 

1491 

0$="  ■ 

1500 

INPUT  "PLOT  KMAX7"  .0$ 

1580 

IF  (J$="V" 

THEN 

GOTO 

1800 

1570 

IF  0$«"N" 

THEN 

GOTO 

1871 

1580 

GOTO  1491 

I 
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i 


1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1671 

1673 

1674 

1675 

1676 

1677 
1673 
1679 
1681 
1682 
1  68  J 
1684 
1690 
1750 
1760 
1770 
1780 
1730 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
I  900 
1910 
1920 
1930 
1  940 
1950 
I960 
1  970 

1  980 

1990 

1991 
2000 
2060 
2070 
2080 
2090 
2100 
2M0 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2201 
2202 


U1 1  =  Lp 

Din  StrS[30] 

StrS=-STRESS  INTENSITY  FACTOR  OIST." 

X  label  5= "L  " 

Ylabel$=''K. MAX- 

CALL  PltlUit  ,Gau5b(  •  )  .Strl  .X labels  .Y  labels  ) 


Ult=Lp 

OUTPUT  2;"K"i 
0$=”- 

INPUT  "PLOT  KMAX-KMIN  AND  R7"  ,QS 

IF  0$="Y"  THEN 

CALL  Dplot I Gaussi •  )  .U1 1  1 

GOTO  1684 

END  IF 

IF  0$=-N“  THEN  GOTO  1684 
GOTO  1674 

I 

QS  =  "  " 

INPUT  "CALCULATE  RATE  UALUES’ ' ,QS 
IF  0S="y"  THEN  GOTO  1780 
IF  0$="N"  THEN  GOTO  2410 
GOTO  1684 

INPUT  “KC  UALUEY  (  MNr!'  3/2  >■  .KI 
IF  L£N(K$)  =  0  Then  GOTO  17S0 
ON  ERROR  GOTO  1780 
K2=UAL<K$I»(  !000".51 
OFF  ERROR 
GCLEAR 

OUTPUT  2; "K"  j 

PRINT  "CRACK  LENGTH’  ,  "RATE.  101^/CYCLE  )" 

PRINT  " 

PRINT 

E^OR  Rt=0  TO  Lp 

Srange“(  GaussIRt  ,2  )+6auss(Rt  ,3)  )/(SOR(  1  .E-3)*2  ) 
5diff*(Gauss(Rt  ,2  )-Gauss<Rt  ,3)  >/(  SOR(  1  .  E-3  )*2  > 

CALL  Rate!  Si'ange  ,Sdi  ff  ,6auss(Rt  ,5  )  ,K2  ) 

PRINT  GaussIRt  ,1  )  ,6auss<Rt  ,6  i 
RpKRl  ,I  )  =  Gauss(Rt  ,1  ) 

RptIRt ,2  )=Gau55(fit  ,6  ) 

File(  2  .Rt  ,1  )=Gauss<Rl  ,1  ) 

File(2  ,Rt  .2  )=Gauss(Pt  ,6  > 

NEXT  Rt 

Fi  let  2 .0 .0 )=Lp 

I 

0$= "  ■ 

INPUT  "DO  YOU  UAr.  :  A  COPY"  , OS 
IF  0$*"Y"  THEN  GOTO  2100 
IF  0$=-N“  THEN  GOTO  2290 
GOTO  1991 

I 

PRINTER  IS  701 

I 

PRINT  Names 

PRINT  "CRACK  GROWTH  RATE  PREDICTIONS." 

PRINT  DATES! TIMEOATE  ) 

PRINT  "RADIUS”  ";  Radius  i  "min .  ” 

PRINT  "APPLIED  STRESS  RANGE'  " : Asmax ; " i Asmin : "MPa .  " 
PRINT  "SPECIMEN  THICKNESS-  "iThS:"  nn." 

PRINT  "SPECIMEN  WIDTH-  "iBSi"  nn." 

PRINT  "A/C  RATIO-  "iAocS 
PRINT  "Klc-  "  !K2rNmn"-3/2' 

PRINT 

PRINT 
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2210  I 

2220  PRINT  "CRACK  LENGTH'  .  "RATE  i  nr  -'CYCLE  )  " 

2230  PRINT  ■■  (nn)“ 

2240  PRINT 

2250  FOR  Pri=0  TO  Lp 

22S0  PRINT  GaussIPri.l  Gauss!  P^;,E) 

2270  NEXT  Pri 
2280  PRINTER  IS  I 

2290  I 

2291  Q$='" 

2300  INPUT  "PLOT7",0* 

23G0  IF  Q$="Y"  THEN  GOTO  2590 
2370  IF  0$=-N"  THEN  GOTO  2410 
2380  GOTO  2291 

2390  CALL  Plt( Lp  ,Rpt( •)  ."CRACK  6R0LTH  RATE ",  "L ",  "RATE "  ) 
2400  I 

2410  I 

2411  0$="- 

2420  INPUT  -STORE  DISTRIBUTIONS’-  , Os 
2480  IF  0*=''Y"  THEN  GOTO  2511 
2490  IF  0$=-N"  THEN  GOTO  2791 
2500  GOTO  2411 

2510  I 

2511  01="- 

2520  INPUT  "CREATE  NEU  FILE’"  ,0* 

2580  IF  0$="Y"  THEN  GOTO  2B40 
2590  IF  Q$*"N'-  THEN  GOTO  2730 
2600  GOTO  251  l 
2510  I 

2520  GOTO  2540 

2530  PRINT  TABXY(4 ,10)j*FILENAME  ALREADY  USED.' 

2640  INPUT  -FILE  NAME’"  .NameJ 

2550  IF  LEN(Narse$  1=0  THEN  GOTO  2540 

2660  ON  ERROR  GOTO  2650 

2570  CREATE  BOAT  NaneS  ,5 

2680  OFF  ERROR 

2690  GOTO  2750 

2700  OUTPUT  2!-K-i 

2710  PRINT  TABXY( 4  , 10  )J "FILE  NOT  FOUND.  TRY  AGAIN,  " 

2720  PRINT 

2730  INPUT  -file  NANE’".Nane$ 

2740  ON  ERROR  GOTO  2700 
2750  ASSIGN  @Pth  TO  Nane$ 

2760  OFF  ERROR 

2770  OUTPUT  @Pth;Pi le( • ) 

2780  ASSIGN  9Pth  TO  • 

2790  I 

2791  0$=-- 

2800  INPUT  "ANOTHER  RUN.  SAME  DIST.-’-.OJ 
2850  IF  OS«"Y"  THEN  GOTO  810 
2870  IF  0$=”N"  THEN  GC^O  2891 
2880  GOTO  2791 

2890  > 

2891  0*="" 

2900  INPUT  "PLOT  GROUP  OP  DISC'S’". OJ 
2950  IF  Q$="Y"  THEN  GOTO  3000 
2970  IF  0$="N-  THEN  GOTO  3020 
2980  GOTO  2891 
2990  I 

5000  CALL  Mplot 
3010  I 

3020  OUTPUT  2:  "K"  i 
3030  GCLEAR 
3040  I 
3050  END 
3050  ' 
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3070  . . . . . 

3080 
3080  I 

3100  SUB  Arrays 
3110  I 

3120  ITHIS  SUBROUTINE  SETS  UP  ARRAYS  UHICH  HOLD  ZEROS  FOR  GAUSS.  INTEG. 

3130  'AND  COEFFS  FOR  C( ALPHA  .BETA > 

3140  COM  /Const/  C<5,4) 

3150  COM  /Zer/  Ero  I  <  S  , Z  )  , Er o2  (  S  .  2  .) 

31G0  RESTORE  3250 
3170  FOR  Uei=1  TO  8 
3180  READ  Ue 
31 80  Erol <  Uei  .  I  '  =  2*Ue 
3200  NEXT  Uei 
3210  FOR  Ro=1  TO  3 
3220  READ  Er 
3230  EroMRo  .2  )=.  '-E^'  2  ■ 

3240  NEXT  Ro 

3250  DATA  0.1894  5'  .0.1  S2S02  .0 . 1  691  ,0 .  149596  .0 . 1  24B29  ,0. 095 1  59  ,0.062254  ,0.02'’' 

52 

5260  DATA  0.095012  ,C. 231604  .0 . 4580 1 7 ,0 . 6 1 7676 ,0.755404  ,0.865631  ,0.54  4  575  ,0.9894 
01 

3270  RESTORE  3360 
3230  FOR  wei=l  TC  6 
3290  READ  Lie 
3300  Ero2(i«ie; 

3310  NEX  *  Ue . 

3320  for  Rorl  TZ  : 

3330  READ  Er 

3340  Ero2-.  Ro  ,  2  r  .  .  R*  ,  5 
3550  NEXT  Ro 

3360  DATA  0. 10122  e.:::36 ,0.31370.0.36268,0,56268 .0.51370  ,0.22238  .0.10122 
3370  DATA  -e , gsCTr  .-0.^9666 ,-0.52553,-0. 18343 ,0. '8343 .0.52553  ,0.79666  ,0.96026 
3380  RESTORE  3440 
3390  FOR  Bb*0  TC  : 

3400  FOR  Aa*i  T2  5 

34  I  0  READ  C(  Aa .Bo  ' 

3420  nE'T  4a 

3430  NE*'  Be 

3440  DATA  0.6164  ,0.04  92  ,-0.4831  .-0. 1746  ,0.7952 .-4.5911  ,17.3181  ,-5C . 5563  ,25 . 287 
7  .-8.4570 ,7.5055  ,-36 .9465 ,75.4833,-75.1 167,26.8666 .-3.8529 
3450  data  2e,E75:  ,-44 , 5540 .44 .2607 ,- 16 .7296 
3460  ' 

3470  SUBEND 

3480  ' 

3490  . . . . . . 

3500  . . . . 

3510  ' 

3520  ' 

3530  SUE  Green :  Grppo  £  _Gf~n  _Lr  1 
3540  ’ 

3550  'THIS  subprogram  EVALUATES  THE  VALUE  OF  THE  GREEN’S  FUNCTION 
3560  I 

3570  'CALC.  C< AlRHA  .BETA  1 
3580  ' 

3590  COM  /Const/  Cl •  1 

3B00  C8b*0 

3610  Alph=1/!1+Lr.) 

3620  FOR  M>1  TO  5 

3630  FOR  N-0  TO  3 

3640  CabrCab  +  C(M,N)«AlpM''(M/2)»Grnpos'"(N/2) 

3650  NEXT  N 

3660  NEXT  M 

3670  I 

3680  ' 

3690  'CALC.  FBETA 
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j700  I 

37 10  Y 1 -ornpos 

3720  Fb  =  <  .2945-.  3912*  Y1  "2+  .7G55*Y1  '■4  -  .  9942  •  Y  1 '  G+  .  5094  •  Y  ?  •&  ) 

3730  I 
3740  I 

3750  iEWLUATE  greens  FUNCTION 
37G0  I 

3770  Grn2  =  50R(2*YI  ) 

3780  Grn3=2»<  1+Fd  .i/S0R(  1  +Y1  )-S0R(2*Yl  ) 

3790  Grn=Grn2+Cab*6rn3 

3800  IF  jYLd  THEN  Grn=Srn/ SQR< I -Grnpos  ) 

3810  SUBEND 
3820  I 

3B30  . . . . . . . . . . 

3840  . . . . . . . . . . . . 

3950 
3BG0  I 

3870  SU6  I  ntegra  1  on(  Lr  .Rad  lus  .Integra  1 5  ,  Integra  12  ,1-d  ,Lcl  ! 

3S80  ' 

3390  COM  /Opt-'  Opt!  ,0pt2 
3900  CON  /Zer/  Er o ! ( .  )  .£r o2 ( *  ) 

3910  CON  /Stress/  Reg,Ual5(lC,20,211,R(10l.BQundC!0,21 

3920  CON  /Cor/  Cor r f lag ,Corr f 1 ag 1 

3930  Integra ! 1 =0 

3940  Integral2*0 

3950 

3950  ‘INTEGRATION  LOOP 
3970  I 

3930  ‘^OR  J='  TO  Ld 
3990  6aussl=0 
4000  6auss2=0 
4010  FOR  Lc:-1  TO  0 
4020  Lc=Lc: 

4030  ' 

4040  IF  J=Ld  THEN 
4050  Usignt=£rol!Lc.l 1 
4000  Yii£rc 1  ( Lc  ,2  I 
4070  GOTO  4120 
4080  END  IF 
4090  Weight=Ero2l Lc  .  1  ) 

4100  Yi=£ro:;Lc 1 
4110  I 
4120  y=Yi 
4  130  I 

4140  'OBTAIN  STRESS  AT  CURRENT  POSITION 
4150  IF  Optl=0  OR  Corrfiag=l  THEN 
4160  C  ipr= Y.Lc 1 /Radius 

4  170  CALL  5 1  ressO'.  C  1  pr  ,  S  t  may  .St  rn  I  n  ) 

4180  GOTO  42G0 
4 1 90  END  I F 

4200  IF  Opt  1  =  1  OR  Opt  1=2  THEN 

4210  Strpo5  =  Lcl«lBouind(J-l  .0)/Lcl+<Y*(  Bound  <  J.0)/Lcl-6ound(J“l  ,0)'Lc!  )); 
4220  IF  J=Ld  THEN  5trpo5  =  Lc  1  » (6ound<  J- I  .0  l.'Lc  1 +<  Y.<  1 -Bound!  J -  1  ,0  1 /Lc  1  1  '  ) 
4230  call  S t  res 5 2 !  S  t  rpos  .S t  max  .S trnin  .Radius  .  J  ) 

4240  I 
4250  END  IF 
4260  ' 

4270  'EUALUATE  GREENS  FUNCTION  AT  CURRENT  POSITION 

4280  IF  Optl-0  OR  Corrflagrl  THEN 

4290  Grnpos'Y 

4300  GOTO  4340 

4310  END  IF 

4320  6rnpo5=Bound(J-l  .0)/Lcl+Y«(Bound(J  ,0)/Lct-Bound(J-l  .0)/Lcl  ) 

4330  IF  J=Ld  THEN  6rnpo5=Bound!  J- I  .0  )/Lc  1 +Y«<  I -Bound! I  .0 '/Lc  1  > 

4340  CALL  Green! Grnpos  ,Grn  ,Lr  ) 

4350  I 
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4550  *5)i-'i*ii’in  1  i  u'l'*  ur  iLAio  ou  f  n/i 

4370  IF  J=Ld  then 

4330  0eno(*>=  I 

4390  GOTO  4420 
4400  END  Ir 

4410  Oenoi<i=SQR<i  I-Bounjfj-I  ,0)/Lct  )/<  BountH  J  ,0  i.'Lc  !  -  Bound!  J  - !  ,0.  ■ .  Lc  ' 
4420  6au35l=S3us5l+Strmav.  •Grn»UeigHt  /Denoi'. 

4430  6au552=Gaus£2-*-Strmi  n*Gr  n*Ue  iQht/Denon 
4440  i 
4  450  * 

44S0  NEXT  Lc2 
4470  I 

4480  IF  J=Ld  Then 

44  90  Integral  I=lntegrall+Gau55l  •S0R(  1  -Bound'  j  -  !  ,0  '  .'Lc  i 
4500  Int  egra  1 2=  I  nt  egra  1 2  +  Gau  ss2»  S0R<  I -Bound’,  J  -  1  ,0!.-uci  i 

4510  GOTO  4550 
4520  END  IF 

4550  I  nt  egra  1  1  =  I  nt  egra  1  '  ♦Gaus  s  I  •SQR(  <  Bound;  J  .0  .‘-Bound'.  J  -  i  ,  2  J  :■  :  .  2 

4540  Integral2=integ’"al2  +  6au5s2»SQR(v  Bound'  J  ,0  i-Bound-  J  -  l  ,C  !  ) '  u  '  ‘2 

4550  NEXT  J 

4560  SL'BEND 

45''0  I 

4580  i 

4590  I* . . . 

4500  . . . . . . 

44 ‘0  1 

4  =  20 

4E50  SUB  Si  f  I  '  Radius  , Gauss!  •)  ,Lg  > 

4  =  40  I 

4550  I  THIS  SUEO.ROGRFiH  CALCULATES  SIF’S  FOP  A  RANGE  C-  CRACK  LENGTHS 
4550  I 

4S'’0  DIM  Cia(  50  ' 

4580  COM  ''Stress  ’'  Reg  ,Ua  1  s (  •  )  ,R(  ♦  >  ,Bound(  •  i 

4590  COH  /Opt/  0ptl,0pt2 

4700  COM  /Aps/  Asna."  .Asnin 

4"’10  COM  /Cor/  Corr  f  J  ag  ,Corr  f  1  ag! 

4720  COM  /Fie/  F;le( .  ) 

4730  COM  /Title/  Th  ,e  ,At- ,T; t le$ 

4740  IF  0pt2=l  THEN 

475'0  PRINT 

4~S0  PRINT 

4770  PRINT  -ENTER  '^ANGE  OF  CRACK.' 

4750  INPUT  "START  M=  RANGE.  (nn)-,Str$ 

4790  IF  LEN(Str$)=0  THEN  GOTO  4770 
4900  ON  error  goto  4770 

4810  Strt=UAL( StrS  ’ 

4820  OFF  ERROR 

4330  INPUT  -END  OF  RANGE.  (nnl’.FiJ 
4340  IF  LEN{Fi*)=0  THEN  GOTO  4850 

4850  ON  ERROR  GOTO  483S 

4850  Fin=UAL(Fj$) 

4870  OFF  ERROR 

4380  IF  Strt’Fin  THEN  GOTO  4770 

4890  INPUT  "INCREMENTAL  LENGTH,  (nn>-,lnl 

4900  IF  LEN<In*)=0  THEN  GOTO  4890 

4910  ON  ERROR  GOTO  4890 

4920  Inc“UAL(In$) 

4930  OFF  ERROR 

4940  END  IF 

4950  IF  Opt2=2  THEN 

4960  La-1 

4970  LOOP 

4980  CIS-'- 

4990  INPUT  "CRACK  LENGTH?  (nn)",Cl$ 

5000  EXIT  IF  LEN(C1$)=0 

5010  Cia'La  )-UAL'Cl$  ' 


TM  MS  1207 


26 


Appendix 


5C:0  La=La+1 
5030  END  LOOP 
5040  Inc=l 
5050  Strt=0 
5060  Fin=La-l 
5070  END  IF 
5080  Corr=0 
5090  OUTPUT  2;  "K"; 

5100  PRINT 
5110  PRINT 

5120  PRINT  "UHICH  CORRECTION  ROUTINE'?- 
5130  PRINT 
5140  PRINT 

5150  PRINT  '■  0.  NO  CORRECTION.  " 

5160  PRINT 

5170  PRINT  ■'  1.  N.  E  R.  WITHOUT  62." 

5180  PRINT 

5190  PRINT  "  2.  NER/SHRIU  HS  CORR." 

5191  0$  =  ''“ 

5200  INPUT  0$ 

5250  ON  ERROR  GOTO  5191 
5250  Corr=UAHOS) 

5270  OFF  ERROR 

5280  IF  Corn',  >0  ftND  Corr?  1  HND  Corr<,2  THEN 
5290  GOTO  5191 

5300  END  IF 

5310  OUTPUT 

5320  IF  Corr«=l  OR  Corr  =  2  THEN 

5530  INPUT  "THICKNESS7  <nn)-.Th 

5340  input  "A/C=  7“  ,Aoc 

5350  Corrflagl=0 

5351  0*=-- 

5360  INPUT  -CO  YOU  WANT  WIDTH  C0RRECTI0N7 ■ ,0$ 

5420  IF  0$=-Y"  then  GOTO  5450 

5130  IF  0$=-N"  THEN  GOTO  5510 

•.4  40  OOT-O  5351 

54-32  Corrflag1  =  l 

'460  INPUT  "WIDTHT  (nnl'.es 

.4"'.'’  ;c  lEN(B$>=0  then  GOTO  5460 

5490  ON  ERROR  GOTO  5460 

j4-i0  B=UAL(B$) 

3=02  OFF  ERROR 
=510  end  if 
5520  I 

5530  OUT^i.T  2!"K-; 

5540  Lp-0 
5550  PRINT 

5560  PRINT  "CRACK  LENGTH"  ,  "KHAX  ,(  MNn'-3/2  )".  "KNIN"  .  "KMAX-KMIN' , ''KrilN/KMAX  " 
5570  PRINT  "  I mn  ) ’ 

5580  ' 

5590  FOR  Lcd  =  ( 5trt  +  Inc  )  TO  Fin  STEP  Inc 
5600  IF  0pt2=l  THEN  Lcl=Lcp 
5610  IF  Opt2»2  THEN  Lci=Cla(Lcp) 

5620  Lc I “GROUND ( Lc I ,81 

5630  Lp=Lp+t 

5640  IF  Opt1“0  THEN 

5650  Ld“l 

5660  GOTO  6040 

5670  END  IF 

5680  FOR  Ld“l  TO  Reg 

5690  IF  Lc  1  <»Bound<  Ld  ,0  )  AND  Lc  I  >Bpund<  Ld- 1  0)  THEN 
5700  SOTO  5730 
5710  END  IF 
5720  NEXT  Ld 

5730  . . TRANSITION  ZONE** . . . 

5740  IF  Lc  I  <Bound(  Ld-1  .0  l-xBoundl  Ld- 1  ,0  )*  I  .5/ 100  THEN 
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B75» 
5760 
5770 
5780 
5790 
5800 
58’0 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
5060 
6070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210 
6220 
6230 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 
6360 
6370 
63o;- 
63- . 
64t' 


Llv  I  j=c:ouna^LU~' 

Ll<2 I'Lcl+Lc I • 1 .5/ 100 
Ll(  3 )=Lcl +Lcl *2/ 100 
FOR  Bioop=l  TO  3 
Lcc=Ll ( Bloop  ) 

IF  B1ood=I  then 

Le=Ld-l 

ELSE 

LB=Ld 

END  IF 

Lr=Lcc/Radiu5 
Corrf lag=0 

CALL  Inlegration(Lr .Radius  ,  integral  I  .  Int egralC  .Le  .Lee  ) 

CALL  Norn! Lee  .Rg 1  .Rg2 , 1 ntegral 1  .Integral2) 

Rg(  I  .Bloop  )=Rg I 
Rg<  2  .Bloop  )=Rg2 
NEXT  Bloop 
FOR  Cod=l  TO  2 

Ddl  =  (Rg(Cod.2)-Rg(Cod.l  l)/(Li;2>-Ll  I  t  i 
Dd2  =  (Rg(Cod.3  )-Rg(Cod  .2  )  )/ (  L 1  <  3  1-Ll  <  2  )  > 

Dd3=(Dd2-Ddl  1/ ( L 1 ( 3  )-L 1 ( 1  )) 

P(  Cod  >=Rg(  Cod.l  )  +  Ddl  ‘(Lcl-Ll!  1  /  '  +  Dc3*'Lcl-Ll(  l  1  )»(Lcl-Ll!2  >  > 
NEXT  Cod 
Lr=Lc 1 /Radius 
Rg1=P( 1  ) 

Rg2=P( 2  ) 

GOTO  6230 
END  IF 


Lr=Lc 1 /Radius 
'CALCULATE  INTEGRAL 
Corrf lag=0 

CALL  Integrati on  (Lr. Radius. Integrall  . Integra  12  .Ld  ,Lc 1  i 

I 

'CALCULATE  NORMALISED  UALUES 

CALL  Norn<  Lr.Rgl  .Rg2, Integrall  .InteQral2) 

I 

'CORRECTION  ROUTINES 

IF  Corr=0  THEN  Cfl=1 

IF  Corr=1  AND  Aoo'=l  THEN 

CALL  Crtnl(Lc1  .Cfl  . T h  .Radi  us ,8 . Aoc ) 

END  IF 

IF  Corr»1  AND  Aoc > 1  THEN 

CALL  Crtn3(Lcl  .Ctl  .Th. Radius, 6. Aoc) 

END  IF 

IF  Corr=2  AND  Aoc<=l  THEN 

CALL  Crtn2(Lc1  ,Cfl  ,Th  .Radius ,8 , Aoc ) 

END  IF 

IF  Corr  =  2  AND  Aoc,-' 1  THEN 

CALL  Crtn4(Lcl  ,Cf  1  ,Th, Radius  ,8  ,Aoc  ) 

END  IF 

I 

Rgl^RgUCfl 
Rg2=Rg2»Cf I 

I 

'RECORD  RESULTS 
Lg=Lp-l 

Gauss!  Lg.l  >“DR0UN0(Lc1  .3) 

Gauss! Lg  ,2  )=0R0UN0! Rg I •Asnax*! PI ‘Lc I • I .E-3 )* .5 ,3  I 
Gauss!  Lg ,3 )=DR0UND! Rg2»A5nax« I  PI ‘Lc I •! .E-3 )' .5 ,3  ) 

Gauss!  Lg ,4 )=Gauss! Lg  ,2  )-Gauss! Lg  .3  ) 

Gauss!Lg  ,5  >“DR0UND!Rg2/Rgl  ,3) 

File!  1  .Lg.l  )“Gauss!Lg.l  ) 

File!  I  ,Lg  .2  I'Gauss!  Lg  ,2  ) 

1 

PRINT  •  "  ;Gauss!  Lg  ,  I  )  ,  ’ "  .Gauss! Lg  .2)  .Gauss!  Lg.5  )  .Gauss!  Lg  ,4  1  .Gauss!  Lg 
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.5  ) 

G4t0  next  Lcp 

G4:0  Fi  le£  1  .0  ,0  )=Lcp-2 

6430  I 

6440  I 

6450  SUBEND 

6460  ! 

6470  I . . . 

6480  I . . . . . 

6490  I 
6500  I 

6510  SUB  StressBl  Strpos  .Strma)^  .Strnin  .Radius  ,J  ) 
6520  Cip=Strpo5 

6530  f  Cl^LC.  STRESS  ACCORDING  TO  INTERPOLATION 
6540  I 

6550  COM  /Stress/  Reg  ,UaI  s  (*)  ,Ri  •)  .Sound*  •  ) 

6560  COM  /Aps/  Asna.H.Asnin 
6570  Y=J-1 

6580  Interp  =  Ual s( Y . I  .R( Y  )  ) 

6590  FOR  X  =  1  TO  ( R( Y  )-I  ) 

6600  E=R(y)-X 

6610  Interp  =  Interp«(Cip-Uai5(Y.E.0'''  +  Ual5iV.I  .E) 

6620  NEXT  X 
6630  Es=Intero 

6640  IF  St rpos I'Bound*  Reg  ,0  1  THEN  E5=0 
6650  ' 

6660  CiDr=Cip/ Radius 

6670  I 

6680  'CALC.  APPlIEO  STRESS 
6690  I 

6700  Ap5=  .5*(  2  +  (  I  /<  l+Cior  )  r2+3«(  i  '  i  I +r  ip-  )  )'4  ) 
6710  I 

6720  :  t'-nax'Aps+Es 

6730  St ’'HI  n“(  Acs*Asnin/Asna.  •‘Es 

6740 

6750  SUBENO 
6760  I 

6770  I . . . . . 

6780  I . . 

6790  I 

6800  SUB  NorntLcI  ,Rgl  .Rg2  , ' ntegra! 1  .Integral2> 
6810  * 

6820  COM  /Aps/  Asnar . Asn 1 n 

6830  Rgl =IntegraI I /PI 

6840  Rg2=Integral2/PI 

6850  I 

6860  I 

6870  SUBENO 

6880  ' 

6890  I . . . . 

6900  . . . 

6910  I 
6920  I 

6930  SUB  StressO* Cipr .Strnax .Strnin > 

6940  I 

6950  COM  /Aps/  Asna.  .Asrtin 

6960  Aps».5*(2+(  I  /  (  l+Cipr  )  r2  +  3.(  l/(  l+Cipr  )  )“4  ) 

6970  I 

6980  Strriax”Aps 

6990  Stmn^Aps.AsMin/Asnax 

7000  I 

7010  SUBENO 

7020  ' 

7030  I . . . 

7040  . . . . . 

7050  I 
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?06W 
7070 
7080 
7090 
7100 
7110 
7!  20 
7130 
7140 
7150 
71B0 
7170 
7180 
7190 
7200 
7210 
7220 
7230 
7240 
7250 
72S0 
7270 
7280 
7230 
7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
7470 
7480 
7490 
7500 
7510 
7520 
7530 
7540 
7550 
7560 
7570 
7580 
7590 
7600 
7610 
7620 
7630 
7640 
7650 
7660 

7670 

7671 
7680 
7740 

7750 


SUB  Plt(Ult.Z:(*>.Str$,Xiabel5  .YlaoelJ  > 

GINIT 

GCLEAR 

GRAPHICS  ON 

Xna;*.  =  Z  2  (1.1) 

Xmi n=Z  2 (  1  ,1  ) 

Vnax  =  Z2  (1,2) 

Ymi n=Z 2 (  1  .2  ) 

FOR  Hh=0  TO  Ult 


IF 

Z2(Hh,1 

)>XWdA 

THEN 

Xnax’sZ z  (  Hh  ,  1 

IF 

Z2(Hh.I 

)  -  Xm in 

THEN 

Xpii  (  Hh  ,  1 

IF 

Z2(Hh  ,2 

) Yf*ia> 

THEN 

YnaA  =  Zz  C  Hh 

IF 

Z2<  Hh  ,2 

)  X  V  w  1  n 

THEN 

Yr»i  n=Zz‘.  Hh  ,2 

NEXT  Hh 

PRINT 

PRINT 

PRINT  -X.RANGE^  ":Xnin:“  TO  ■■:Xnax 
PRINT  "Y. RANGE-  ''sYnin:"  TO  ■■:»nax 
INPUT  "X-AXIS  label  SPACING’"  .XapI 
IF  L£NfXsp$)-0  THEN  GOTO  7250 
ON  ERROR  GOTO  7250 
X5BC=VAL( XscS  ) 

OFF  ERROR 

INPUT  "Y-AXIS  LABEL  SPACING’’  .Vso* 

IF  L£N(Xsp*>=0  THEN  GOTO  7300 
ON  ERROR  GOTO  7300 
Yspc=UAL( YspS  ) 

OFF  ERROR 
OUTPUT  2!"K"i 
MOUE  35,95 
LABEL  Str$ 

UIEWPORT  15,115,30,90 
FRAME 

Xdiff-Xrnax-Xmn 
Ydi f  f -Ymax- Yni n 

WINDOW  Xnin.X«ia.<  +  Xdiff/t8,Yrnn-Ydiff/l0,YnaA  +  Ydiff/l0 
AXES  X ape /5 , Yspc /S  ,0  ,0  ,5  ,5 , 1 
CLIP  OFF 
C51ZE  3 

FOR  Xlab=0  TO  Xmax tXdi f f / 1 2  STEP  Xspc 
IF  XlabOimin-Xdif  f/12  THEN  GOTO  7500 
MOUE  Xlab-Xdiff/140 , Yn i n-2» Ydi f f / 1 0 
LABEL  USING  "#  ,K- ;Xlab 
NEXT  Xlab 

MOUE  Xi*tax  +  >  di  f  f /8  ,Ynin-Ydi  f  f /S 
LABEL  USING  ' «  .K  "  : X labe 1$ 

FOR  Ylab=0  TO  Ynax  +  Yd i  f f / 1 2  STEP  Yspc 
IF  YlabXYnin-Ydif f/12  THEN  GOTO  7570 
MOUE  Xnin-Xdif f/10  .Ylab-Yd:f f/30 
LABEL  USING  "»  ,K" ;Ylab 
NEXT  Ylab 

FOR  Ylab=0  TO  Ynin-Ydi f f / 1 0  STEP  -Yspc 
MOUE  Xnin-Xdiff/l0,Ylab-Ydiff/30 
LABEL  USING  "#  ,K"  ;Ylab 
NEXT  Ylab 

MOUE  Xnin-Xdi f f / 10  ,Ynax+Ydi f f /8 
LABEL  USING  " J  ,K " ; Y labe  1  $ 

P-1 

FOR  Rr-0  TO  UH 

PLOT  Z2(Rr,l  ),Z2(Rr,2),P 

NEXT  Rr 

0*-'" 

INPUT  •COPY’",0* 

IF  0$-'Y"  THEN  GOTO  7780 
IF  0$-"N"  THEN  GOTO  7800 
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77B0  GOiO  7B71 
7770  I 
7780  OUTPUT 
7790  I 
7800  SUBEND 
7810  I 

7820  I . . . 

7850  I . 

7340  I 

7850  SUB  I npu t 5 ( Rad lus ) 

7860  I 

7870  COM  /Opt/  Opti .Opt2 

7880  COM  /Stress/  Reg  .Ua 1 s < » >  ,R< ‘  i  ,8ound< •  ) 

7990  COM  /Aps/  Asrax.Asnin 
7900  COM  /Fie/  Fi  lef •  ) 

7910  COM  /Title/  Th  ,B  .  Aoc  ,Nane$ 

7920  DIM  Lr(50.2  i  .!:<  102.2) 

7930  PRINT 
7940  I 

7950  IF  Opt  1=2  THEN 

7960  PRINT  "RSF  LOADED  FROM  FILE- 

7970  GOTO  8000 

7980  PRINT  "FILE  NOT  FOUND." 

7990  OFF  ERROR 

8000  INPUT  "NAME  OF  FILE7-  ,Name$ 

8010  IF  LENl Names  )=0  THEN  GOTO  8000 
8020  OUTPUT  2i"K"i 
8030  PRINT 
8040  PRINT 

8050  ON  ERROR  GOTO  7980 
60B0  ASSIGN  @Path  TO  Names 
8070  OFF  ERROR 
8080  ON  ERROR  GOTO  8100 
8090  ENTER  @Path;lz( •  ) 

8100  OFF  ERROR 

8110  ASSIGN  @Path  TO  • 

8120  PRINT  "FILE:  "  .NameS 

8130  PRINT  “NO.  OF  ENTRIES:  "  ,L3<2  ,0  ) 

8140  PRINT  "RADIUS  OF  HOLE : “  ,L: ( 3  ,0  ) 

8150  Radius=Lz< 3  ,0  ) 

8160  UIt=L2( 2  .0 >-l 
8170  I 
8180  I 
8190  ELSE 

8200  Name$="KEyEOARO  ENTRY." 

8210  OUTPUT  2:  "K"  i 
8220  INPUT  ■■RADIUS')  (nnl'^.RadiS 
8230  TF  LEN(RadiS)=0  THEN  GOTO  8220 
8240  ON  ERROR  GOTO  8220 
8250  Radius  =  UAL( RadiS  ) 

8260  OFF  ERROR 

8261  PRINT  "ORIGIN  OF  DISPLACEMENT  COORDS." 

8263  0$="" 

8270  input  "0  FOR  CENTRE  OF  HOLE.  1  FOR  E0GE.“.0$ 
8330  ON  ERROR  GOTO  8263 
8340  Kd=VAL(0$> 

8350  OFF  ERROR 

8360  IF  Kd<>0  AND  Kd/M  THEN  GOTO  8263 
8370  Lr(l,0)"Kd 
8380  I 

8390  L;(  3 ,0  )*Radius 

8400  PRINT  "INPUT  POSITIONS  AND  STRESS  VALUES.” 
8410  PRINT 
8420  PRINT 

8430  PRINT  "NO."  ."X."  ."STRESS." 

8440  Inp=-I 
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845a  Luur 
8460  Inp=Inp+l 
8470  Ual;'.$="" 

8480  INPUT  "UALUE  OF  X  T  (r>m>“,Ualx$ 

8490  EXIT  IF  LEN<OalxI)=0 
8500  ON  ERROR  GOTO  8480 
8510  L;(  Inp,!  )=UAL<Ualx$) 

3520  OFF  ERROR 

8530  INPUT  “UALUE  OF  STRESS’  (MPa)",Uals$ 

8540  IF  LEN(Ual5S)=0  THEN  GOTO  8530 
8550  ON  ERROR  GOTO  8530 
85G0  Ualstr=yAL< UalsS ) 

8570  OFF  ERROR 

8580  Lr( Inp  .2  )=Ualstr 

8530  PRINT  Inp  ,L:< Inp  .  I  )  .L'( Inp  ,2  ) 

8600  END  LOOP 

8601  0$='"' 

8610  INPUT  "CHANGE  ANY’‘.0$ 

8670  IF  0$="Y"  THEN  GOTO  8710 
8680  IF  Q$="N"  THEN  GOTO  87G0 
8690  GOTO  8601 

8700  I 

8710  INPUT  -ENTRY  TO  BE  CHANGED’”  .Nn 

8720  INPUT  “UALUE  OF  X.  ( mn  )  ”  .Lz< Nn  ,  1) 

8730  INPUT  -UALUE  OF  STRESS.  ( MPa  )" ,Lz( Mm  ,2 > 

8740  PRINT  TABXYC  I  .Nn  +  G  1  ,Nn  ,Lz( Nn  ,  1  ) ,Lz(  Nn ,2  ) 

8750  GOTO  8601 

8760  Ult=Inp-l 

8770  Lz! 2 .0 )«Ino 

8780  1 

8781  0$'”" 

8790  INPUT  -STORE  DISTRIBUTION’"  .0$ 

8850  IF  0S=-Y-  THEN  GOTO  8881 
8860  IF  0$=-N-  THEN  GOTO  9190 
8870  GOTO  8781 
8880  1 

8881  0S=" 

8890  INPUT  -CREATE  NEW  FILE’-.0$ 

8950  IF  0$--Y-  THEN  GOTO  8990 
8960  IF  0$--N-  THEN  GOTO  9100 
8970  GOTO  8900 

8980  1 

8990  GOTO  9010 

9000  PRINT  TABXYI4.10);-FILEN- ME  ALREADY  USED.  TRY  ANOTHER." 

9010  INPUT  "FILE  NAME’-. Names 

9020  IF  LEN(Name$)=0  THEN  GOTO  9010 

9030  ON  ERROR  GOTO  9000 

9040  CREATE  BOAT  NameS  ,5 

9050  OFF  ERROR 

9060  GOTO  9120 

9070  OUTPUT  2;-K-i 

9080  PRINT  TABXY(4,10)i-FILE  NOT  FOUND.  TRY  AGAIN." 

9090  PRINT 

9100  INPUT  "FILE  NAME’", NameS 
91.0  ON  ERROR  GOTO  9080 
9120  ASSIGN  @Pth  TO  NameS 
9130  OFF  ERROR 

9140  OUTPUT  ePthjLzI •  ) 

9150  ASSIGN  gPth  TO  • 

91G0  ' 

9170  END  IF 

9180  GOTO  9234 

9190  I 

9200  INPUT  "TITLE’"  .Names 
3201  0$-"" 

9210  INPUT  -COPY  OF  COORDS’"  .0$ 
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a::n  if  qs^-y”  iHtN 

92IZ  PRINTER  IS  70! 

9215  PRINT  NaneS 

3214  PRINT  "RflOIUS'"  jRadius  : nn.‘ 

9215  IF  Lz{  I  .0)=0  then 

9216  PRINT  "L  MEASURED  FROM  CENTRE  OF  HOLE.’ 

9217  END  IF 

9218  IF  Lz(  I  ,0)=l  THEN 

9219  PRINT  "L  MEASURED  FROM  EDGE  OF  HOLE." 

9220  END  IF 

9222  PRINT 

9223  PRINT  "L  <  mn  STRESS  (MPa)" 

9224  FOR  Po=0  TO  Ult 

9225  PRINT  L2< Po  ,  1)  .L;( Po  .2  ) 

9226  NEXT  Po 

3227  PRINTER  IS  1 

9223  GOTO  9234 

9230  END  IF 

9231  II  THEN  GOTO  9234 

3232  GOTO  9201 

9235  I 

9234  IF  L;<  1  ,0  )=0  THEN 

9240  FOR  Yy=0  TO  50 

9250  L;( Yy  ,  1  )  =  L;( Yy  .  1  (-Radius 

9260  NEXT  Yy 

3270  END  IF 

9230  call  PIKUlt  ,Lz(- ).NaneI."X"  ,''STR"  ) 

5290  INPUT  "’NO.  OF  REGIONS’ "  ,Re$ 

9300  IF  LEN<Re$)=0  THEN  GOTO  923C 

9310  ON  ERROR  GOTO  3290 

9320  Reg=UAL(Re$) 

9330  OFF  ERROR 

9340  GCLEAR 

9350  FOR  Z=0  TO  Uit 

5350  PRINT  Z  ,L2(  Z  .1  )  ,Lz(  Z  ,2  ) 

9370  FiIei0,Z,l  )=Lz(Z,t  ) 

9390  File(0,Z  .2  )=Lz(  Z  .2  ) 

9390  NEXT  Z 

3400  File(0 ,0  .0  )=Ult 

94 1 0  8oun( 0  )=0 

9411  OS=’'' 

9413  INPUT  "COPY.-, OS 

9414  IF  0J="Y''  THEN 

9415  PRINTER  IS  701 

9416  PRINT  Names 

9417  PRINT 

9419  PRINT  -L  (mm)"  ."STRESS  (MPa)" 

9420  FOR  Pp=0  TO  Ult 

9421  PRINT  Lz(Pp  ,1  )  ,Lz(Pp  ,2  > 

3422  NEXT  Pp 

9423  PRINTER  IS  I 

9424  GOTO  9428 

9425  END  IF 

9425  IF  OS-"N"  THEN  GOTO  9428 

9427  GOTO  941 1 

9428  PRINT  ’INPUT  NUMBERS  OF  POSITIONS  OF  THE  BOUNDARIES." 

9430  PRINT 

9440  PRINT 

9450  PRINT 

9450  FOR  X=^I  TO  Reg 

9470  print  USI’G  ,K’ ;  "POSITION  OF  ENO  OF  REGION’./ 

9480  INPUT  Bour ( X  ) 

9490  Bound! X  ,0  )“Lz( Boun( X  )  ,  I  ) 

9500  NEXT  X 

9510  I 

9520  Bound!  0  ,0  )=Lz{0  .  I  ) 
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9b5W  fUK  onos=«;  lu  ncy-i 
3540  Lpn=0 

9550  FOR  Pt s=Boun( Bnds )  TO  Boun<Bnd5+I' 

95G0  Lpn=Lpn+) 

9570  Oals( Bnds  .Lpr  ,0  )=L2(Pt5  .1  ) 

9580  0al5<Bnd5,Lpn,l  )-Lz(Pts,Z)/Aanax 

9590  NEXT  Pts 

9G00  R(Bnd5)  =  Lpn 

9B10  NEXT  Bnds 

9BZO  ! 

9630  'CALC.  DIUIDED  DIFF.  TABLE 
9G40  FOR  Z1=0  TO  (Reg- I  ) 

9G50  FOR  Xl=2  TO  <R(  Zl  I-*-!  ) 

9BG0  FOR  Yl  =  l  TO  ( R< Z I  )-( X  I  -  1  ) ) 

9G70  Ii-(Uais(ZI  ,YI  +  I  ,X1-1  )-Ual5(Zl  .YI  .Xl-I  ))/(0al5(ZI  ,YI+X!-I  ,0)-Oals(2l  ,v|  ,i2 
)  ) 

9G80  Oaj5(  Z  I  ,V  1  .X  n-DROUND(  1 1  .5  ) 

9B90  NEXT  Yl 
9700  NEXT  XI 
9710  NEXT  Zl 
9720  I 

9730  'PLOT  routine 

9740  I 

9750  Ult=0 

97BO  yy=0 

9770  Xni n=Bound ■ 0  . 0  I 

9780  Xnax=Boura  >-:  .0) 

5790  X d 1 T f  =  X na » - > n  1  n 

3800  FOR  Range’Xnin  TO  Xna*  STEP  (XdiffT|00) 

9810  Ult»UIt+l 

9820  IF  Range 'Bound' Yy 1 1  ,0  )  THEN  Yy  =  Yy+I 
9630  Interp  =  uai»i'y,l  ,R(Yy)) 

9840  I 

9850  FOP  Xx. I  TC  • F! Yy  > 

98B0  E  =  R( Yy  i-x  . 

9870  I  nt  er-px  I  n  t  erD»  (  Range-Ual  s{  Yy  .E  .0  )  )+Vals(  Yy  ,  1  ,E  ) 

9880  NEXT  Xa 

9890  Zr( U1 1  , 1  ixPange 

9900  Zz(  U1 1  , Z  1=  1  n t  erp ♦  Asmax 

9910  NEXT  Range 

9920  Din  Strstsei 

9930  Str$x'RESIDUAL  STRESS  FIELD' 

9940  CALL  PI t ( U1 1  .Zr( •  )  ,StrS  .'L' ,"STR‘ ) 

9950  I 

9951  0$-" 

9960  INPUT  ■CK-’-  ,0S 

10020  IF  OSx-Y-  Then  60TO  10090 

10030  IF  OSx-N"  THEN 

10040  U1  t=Lz(  2  ,0  >- I 

10050  GOTO  9280 

10060  END  IF 

10070  GOTO  9951 

10080  I 

10090  SUBEND 

10100  I 

10110  . . . . . . . . 

10120  I . . . . . . . . . 

10130  I 

10140  SUB  Rate(U,W,R2  ,K2  ) 


. . . . . . 

101B0I . •••♦CRACK  RATE  DATA  :  NAME  ,KCRIT  .DATA 

|0170I^^^^^^^*^^^^»^^^^^^(ALL  in  . . 

101801 . •••• . . . 


10190  DATA  L71 
10200  KcriLd«l .86E+3 
10210  Kl-Kcntd 
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I022B  -2,yB.a,l02,llb.lA/.jlb  .bbfc' ,  I  BUU  ,-l  .It-/-  ,i  .b£-o  ,4t- 

E-2,0 


10230  DATA  -I  ,65.9  ,69.0  .80.0 , 101  .455 ,1  .IS  .-!  .tE-7  .4E-7  .lE-S  .36-6  ,4E- 
10240  DATA  -.5  ,49.5.51  .0  .65.0 .450.-’50.-l  .  IE-7  .4E-7 .1  .5E-6  ,  IE-3  ,  IE-2  , 
10250  DATA  -.25  .41 .3 ,43.0  ,50.0  .450  ,620.-1  .IE-7 .4E-7  .9E-7,2E-3  ,  lE-2  ,0 
10260  DATA  0 ,33.0 ,55.0  ,40.0  ,390 ,510  I  .IE-7  .3E-7  ,7.3E-7  ,2E-3  ,  i'e-2  ,0 
10270  DATA  .25 ,26.6  ,29.0  .52.0  ,235  ,390,-1  .IE-7  .2E-7  ,4E-7  ’6E-4  ,'lE-2  ’o 
10280  DATA  .5,23.4.24.5,28.0.149,300.-1 , IE-7. 2E-7  4E-7  1.55E-4  IE-' 
10290  DATA  999 

. . . . . . . . 

10310  Opdev=701 
10320  Data$="YES" 

10330 

103401 . ••‘•‘READING  IN  OE  . . . 

103501  . . . . . . . . 


. . . 

10370i““CRACK  RATE  DATA“** 

103801  ••* . . . 

10390  READ  RatenaneS 
10400  11=0 
10410  LOOP 
10420  READ  Ts 
10430  EXIT  IF  Ts=999 
10440  R<  1 1  )=T5 
10450  12=0 
I04B0  LOOP 
10470  READ  Dl; 

10480  EXIT  IF  Dk=-1 .E+0 
10490  Ak  (  II  ,12  )=L0G(DI'  > 

10500  I2“I2+I 

10510  END  LOOP 

10520  13=0 

10530  LOOP 

10540  READ  Rte 

10550  EXIT  IF  Rte=0 

10560  6(11  ,15  )=LOG(Rte ) 

10570  13=13+1 
10580  END  LOOP 
10590  11=11+1 
10600  END  LOOP 
10610  1=11-1 
10620! 

10630  IF  U+U>K2  THEN  GOTO  10850 
10640  IF  U^=0  THEN  GOTO  10830 
10650  IF  u+U(=0  THEN  GOTO  10850 
10660  R3=(  U-Ui  )/(  U+W  ) 

10670  P=U+W 

10680  IF  R3;'R(0)  then  GOTO  10870 
10690  R3=R(0) 

10700  Ul •LOG<P/2«( I-R3  )  1 
10710  J=0 

10720  IF  Ul<AE(J,0i  THEN  GOTO  10830 

10730  FOR  K-l  TO  9 

10740  IF  Ak(J,K)^WI  THEN  GOTO  10770 

10750  IF  Ak(J.K+n«0  THEN  GOTO  10770 

10760  NEXT  K 

10770  KO'UI-AEI J  .K  ) 

10780  Z-S(  J  ,K  )+K3‘(S(  J  .K  >-S(  J  .K-l  >>/<  AKt  J  .K  >-Ak(  J  .K-1  )  > 
10790  IF  K3<0  THEN  GOTO  10810 

10800  Z-Z+K3»K3/( (LOG( K1 •( I-R5  1/2  )-U! +K3 ) ■2-K3»K3 ) 

10810  R2=EXP(  Zl^lll-P/Kl  )/<  I-<U+U  )/K2  ) )’  .5 

10820  GOTO  I  1250 

10830  R2=0 

10840  GOTO  11250 

10850  R2=l .E+10 

10860  GOTO  1 1250 


,  4  fc  -  b  .  1  ,  I 

4 , IE-:  ,0 
0 

0 
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t08<'W  WI=LUtnw; 

10880  IF  R(I)>R3  THEN  GOTO  10930 
10890  R3=R< I  ) 

10900  J=I 

10910  P=2*W/( 1-R3) 

10920  GOTO  10720 

10930  FOR  J=l  TO  I 

10940  IF  R<J>>R3  THEN  GOTO  109G0 

10950  NEXT  J 

10960  F=<R3-R(  J-l  >  )/(R(  J  )-R(  J-  t  ) ) 

10970  H3=F»flk(  J  ,0  >+<  1-F  >*ftk(  J-]  .0) 

10980  IF  U1<H3  THEN  GOTO  10830 
10990  K=0 
1 1 000  L=0 
1 1010  H4=S( J ,0 ) 

I  1020  IF  S(  J  .K+1  )>S(  J-l  .L+1  )  THEN  GOTO  11140 
I  1  030  K=K+ 1 

11040  IF  S( J  .K  )=S( J-l  ,L  )  THEN  GOTO  11020 
I  1050  Hl=St  J  ,K  ) 

1  I0B0  G=(S(  J-l  ,L+n-S(  J-l  .L  )  )/(  Ak(  J-l  .L+1  )-AK  J-1  ,L  )  ) 

1  1070  H2=Ak<  J-I  .L  )+(S(  J,K  )-S(  J-I  ,L))/G 

11080  H5=F*Ak(  J  .K  )+<  l-F  )»H2 

11090  IF  Ak<J.K+1)=0  THEN  GOTO  11220 

11100  IF  UKHS  THEN  GOTO  11220 

11110  H3=H5 

11120  H4=HI 

11130  GOTO  1 1020 

11140  L=L+! 

11150  IF  S( J  .K  )=S( J-l  .L  )  THEN  GOTO  11020 
11160  HI =S( J-l  ,L  ) 

11170  G=(  S(  J  ,K+1  >-S(  J  .K  )  )/<  Ak<  J  ,K+I  )-Ak(  J  ,K  )  1 
I  I  180  H2=Ak(  J  .K  )+(  S(  J-l  ,L  )-S(  J  ,K  )  )/6 
11190  HS»F*H2  +  (  1-F)*Ak(  J-l  .L) 

11200  IF  Ak( J-l  .L+l  )=0  THEN  GOTO  11220 
11210  GOTO  1 1100 
1 1220  K3»W!-H5 

11230  Z=H1+K3«(HI-H4  )/(H5-H3 ) 

1 1240  goto  10790 
1 1250  R2=DR0UND(R2 ,3) 

1 1260  SUBEND 
I  1270  I 

1  I  280  I  . . ...» . . . . . 

1  1290  I . . . . . . . 

1 1 300  1 

11310  SUB  Crt  n  1  (  0  ,Fch  ,  T  ,Rad i us  ,B  , Aoc  ) 

11320  con  /Cor/  Corr f i ag  ,Corr f 1 ao I 
11330  I 
I  1 340  Thi=0 
I  1 350  A=C»Aoc 

1 1360  0-1  +  1 .464»( A/C  I"! .65 
1 1370  Ml-1 . I3-.09»(A/C  ) 

11380  M2--.54+.e9/<  .2  +  ( A/C  ) ) 

1  1  390  M3-.5-I /(  ,G5+( A/C  )  >+l 4» ( 1 -A/C  1 '24 
I  1400  6l-l  +  <  .  l  +  .35»<  A/T  r2  )•(  l-SIN(Thi  )  )"2 

11410  G3-M  +  .04»<  A/C  )  )»<  1  +  .  1  ♦<  l-C0S(Thi  )  )'2  >•(  .85+.  15*  (A/T  )*  .25  ) 

1  1420  Fthi-((A/C)''2«C0S<Tln  )-2  +  SIN(Thi  1*2)". 25 
11430  IF  Corrflagl-1  THEN 

11440  Fu-( l/COS(PI*Radius/(2»B  ))»l/C0S(PI*(2*Radiu5+C  )/(4«(B-C )+2*C  )*S0R( A/T  )  )  ) ' 
.5 

1 1450  ELSE 
1 1460  Fw-t 
1 1470  END  IF 

I  1480  Fch-(  (ni+M2*(A/T  )'-2+f13»(  A/T  )'  4  )*G  1 ‘GS-F  t  hi  *Fu>  )/S0R(  0  ) 

11490  K I -S0R( ( 4/PI+A»C/(2*T«Radlus  1  )/( 4/PI+A»C/<T •Radius  )  )  ) 

11500  Fch-Fch«KI 
I  1510  5U8END 
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I  I  520  I 
11530  I 

1  IS40  . . . . . . . . . 

11550  I#**** . . . . . . 

1 1560  I 

11570  SUB  Crtn2< C  .Cf!  ,T  .Radius  ,B,Aoc  ■ 

11580  ' 

1  1 530  I 

1  1  500  COM  /Cor/  Corrf 1 ag  ,Ccrr f  1  ag 1 
I  1610  Thi=0 
1 1 620  A=C*Aoc 

1 1630  0=1  +  1  .464*<  A/C  )' I .65 
1 1640  M1=1 . l3-.09«( A/C ) 

11650  M2“-.54+.89/<  .;  +  (A/C  1) 

I  1660  f13=.5-l/<  .B5  +  (A/C  >  >+ 1  4* <  A/C  )“24 
1  1670  61  =  l+(  .  l  +  .35«vA/T  )“2  )•(  1-SIN(Thi  )  )-'2 
1  1680  Lam=1/(  l+(C/Radias  )*CCS( ,85*Thi )) 

11690  G2=!  1+  .  358«4ar»  ;  .  425* Law -'2-  1  .578* tan' 3+2 . 1 56*Lam"  4  )/(  1  +  .  1  3«Lari  2  ) 

1  1  700  63=  (  l  +  .04  +  A/C‘*'  :  +  .!•'  1  -C0S(  Thi  )  )-2  )•<  .85+,  15»(A/T>".25) 

11710  Fthi  =  ((A/C''2»C0S(Thi  1‘2+SIN(Thi)“2)'  .25 
11720  IF  CorrfiaQl=I  ThEN 

117o0  6w=(l/C0SlFI‘ Radius /(2»B))*l/CGS<PI  +  i2«Racliu5  +  C)/(4*(8-C)  +  2»C)»SQR  (A. 'T))) 

1 1740  else 
1 1 750  Fw=l 
1 1760  END  IF 

I  1  770  Nr  =  <  <  Ml  +M2*  +  '  2»M3*(  A/T  )•  4  ' *5 1  •S2  +  63*F t h i ‘F w  )/SQR(Q  ! 

1  1  780  K  1  =SQR(  (4/4 ;  •M«:/!2«T*Radius))/(4 /PI +A*C/(T  + Radius))) 

1 1 730  Nr=Nr*K 1 
1 1 800  I 
11810  I 

1 1820  Lr=C/Radiu; 

1 1 830  Corr  f !  35=  I 

11840  CALL  Integral . O'M  Lr  .Radius  ,  Integral  I  ,  IntegralZ  .  1  .C  ) 

11850  CALL  Norr  . L+ .4g '  ,Rg2  .  Integral  I  ,Integral2  ) 

1 1 860  Shi  V 1 =Rg  ! 

1  1870  Cf l=Nr/5hivl 
1  1880  I 

11890  SueSflD 

1  1 900  I 

11910  I . . . . . . . . . 

1 1920  I . 

11930  Sue  Cr  t  n3(  C  ,Fc  n  .  T  ,Rad  1  US  ,B  ,Aoc  ) 

11340  COM  /Cor/  Co-rf  lag  .Corrf  I'agl 
1 1 950  I 
1 1960  Thi=0 
1 1 970  A=C«Aoc 

I  1 980  0=1  +  1 . 464. ( C/A  <"1.65 
11990  Ml  =S0R<  C/A  )•!  I  .  .04.<  C/A  )  ) 

I  2000  n2= .2* ( C/A  ;'4 
12010  M3='.  1  l.(C/A  )-4 

12020  Gl  =  l  +  i  .1  +  .35*(C/A).(A/T)"2)»(l-SIN(Thi))"2 

12030  G3  =  (  1  .  13-. 0g»(C 'A  1  ).'  1  +  .  ].(  I  -cost  Thi))'2)*(.85+.15.<A/T)''  25) 

12040  Fthi  =  <  (C/A)-2*SlMThi  )'2  +  C0S(Thi  )"2)'-.25 
12050  IF  Corrflag!=l  THEN 

12060  Fu=(  l/C0S<PI«Radiu5/(2*B))+l  /COSt  P I  •(  2»Radiu5+C  >/<4*(B-C)+2»C  )»S0R(  A/T  > )  )' 

12070  else 
12080  Fu-1 
12090  END  IF 

12100  Fch=( (MI+M2  +  <  A/T  )-2  +  M3*( A/T  >-4  >*6 1  »G3*F t hi •Fu )/SOR< 0 ) 

12110  F I-SORl <4/PI+A»C/<2»T«Radius)>/(  4/PI+A.C/(  T  +  Radius  )  )  ) 

12120  Fch=Fch+Kl 
12130  SUBEND 
12140  I 

12150  I . . . 
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12lbli)  . . . . . . . . 

12170  I 

12180  SUB  Crtn4(C  ,Cf I  .1  .Radius  .8  ,Aoc > 

12190  COM  /Cor/  Corrf lag .Corrf lag! 

12200  Thi-0 
12210  A=C»Aoc 

12220  0=1  +  1  .4S4*(C/A)'  1  .65 
12230  Ml =S0R(C/A )•( I  +  .04  +  (  C/A ) ) 

12240  M2=.2*(C/A )-4 
1 2250  M3=-. 1 1 •(C/A )-4 

1  2260  ei  =  l  +  (  .  I  +  .35»(  C/A  >*(  A/T  >''2  )•(  l-SIN(  Thi  >  r2 
12270  Lar>i=l/!  1+(C/ Radius  )»C0S(  .85*Thi  >) 

12280  62  =  <  1  +  .358*Lan+l  .425»Lai>i“2-l  .578*Lan~3+2.  l56»Lan“4)/'(  1  +  .  13'Lat"  2  ' 

12290  G5=(  1  .  I5-.09*C/A  )•(  1  +  .  1*(  1-C0S(  Thi  )  >“2  >•(  .05+.  I  5«<  A/T  )•  .  25  ' 

1 2300  Fthi  =  ( (C/A  )“2*5IN( Thi  )"2  +  C0S(Thi )“2 )' .25 
12310  IF  Corrflag1=l  THEN 

1 2320  Fu=( 1 /C0S( PI*Radiu5/( 2*8 )  )•! /C0S<PI*(  2»Radius+C  )/( 4*<  6-C  )  +  2*C  )»S0R<  A  T ) . 

.5 

12330  ELSE 
12340  Fu=! 

12350  END  IF 

12360  Nr=( ( Ml +M2»( A/T ) ' 2  +  M3* ( A/ T ) ' 4 ) 'G 1 •G2»G3*F t h i ‘Fu  )/S0R( Q  ) 

1 2370  K  t=S0R( <  4/PI+A»C/( 2*T*Radiu5  )  )/<  4/PI +A»C/ <  T 'Rad lus  )  )  : 

1 2580  Nr=Nr*K I 
12390  I 
12400  I 

12410  Lr=C/Radius 
12420  Corrflag-1 

12430  CALL  Integral ion( Lr  .Radius  .Integral  1  ,Integral2  , 1  ,C  ) 

12440  CALL  Norn! Lr  .Rgl  .Rg2  .  Integral t  ,  Integral2  ) 

12450  Shivl=Rgl 
12460  Cft=Nr/Shivl 
12470  I 
12480  SUBEND 
12490  ! 

12500  I . . . 

12510  . . . . . . .  . 

12520  I 
12530  I 

12540  Sub  Mpiot 
12550  I 

12563  DIM  MplK  5 .2  .50  .2  ) 

12570  DIM  Hid(2  .50  ,2  ) 

1  2580  DIM  Pid(5  .50  .2  ) 

12590  I 

12600  GCLEAR 

12610  OUTPUT  2rK"; 

12620  I 

12630  Nn=0 

12640  I 

12650  LOOP 

12660  Nn*Ni>i+  I 

12670  OUTPUT  2;-K"i 

12680  GOTO  12710 

12690  OUTPUT  2i”K";  J 

12700  PRINT  "FILE  NOT  FOUND." 

12710  INPUT  -NAME  OF  FILE’-.NameS 
12720  IF  LEN(Name$)-0  THEN  GOTO  12710 
12730  ON  ERROR  GOTO  12700 
12740  ASSIGN  ePathI  TO  NaneS 
12750  OFF  ERROR 
1  2760  ENTER  SPathl iHld< •  ) 

12770  ASSIGN  @Pathl  TO  • 

12780  FOR  X»0  TO  2 
12790  FOR  Y«0  TO  50 
12800  FOR  Z«0  TO  2 
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12810  Mpl  t  (Nn  ,X  ,y  ,£  >“HldC  X  ,y  .Z  ) 
1:820  NEXT  : 

12830  NEXT  Y 
12840  NEXT  X 
12850  I 

12860  I 

12861  0S='‘" 

12870  INPUT  "ANOTHER  FIlE’-.QS 
12930  IF  QS="Y"  THEN  GOTO  12970 
■2940  EXIT  IF  OS="N“ 

12350  GOTO  12861 


12960 
1  2970 
‘2S80 
1  2990 

■  3000 
13010 
13020 
:  3030 
1  3031 

■  3840 
13100 
13!  10 


END  LOOP 

PRINT  "C.  EXIT." 

PRINT  "1.  residual  stress  FIELD." 

PRINT  "2.  STRESS  INTENSITY  FACTOR  DIST. 
PRINT  "3.  CRACK  GROWTH  RATE  DIST." 

Nd=t 
0S="  • 

INPUT  "WHICH  DO  YOU  WANT  TO  PLOT’"  ,0S 
ON  ERROR  GOTO  13031 
Nd=UAL( OS  ) 


'5120  OF-  ERROR 

'3138  IF  Nd-  '0  AND  Nd  ,  1  AND  Nd  -2  AND  Nd-  3  THEN 
13140  GOTO  13031 


13153  END  IF 

■3150  IF  Nd=0  THEN  GOTO  13450 
'3170  FOP  Hd=l  TO  Nn 
'3!Se  FOR  Xr=0  YO  50 


:  31  93  Pldi  rid  ,Xr  ,  I  )=Ad1  t  <Md  ,Nd-  I  .Xr  ,  1  ) 

■  3200  Pld'  Md  ,X;  ,2  )  =  NdI  t  f  Md  ,Nd- 1  .X:  ,2  ) 
13210  rid!  Md  .0  ,0  i^MdI  t  (  Md  ,Nd- I  ,0,0) 

' 3220  NEXT  xj 
1  3230  NEXT  fid 
'3240  DIN  LabS[301 
'3250  IP  Nd=I  THErj 
'32G0  XS="X'' 

13270  yS=“OTR" 

•3:38  ■^abS=''RE5!DuAL  STRESS  DIST." 

■  32  90  END  IF 
13300  IF  ^a-C  Then 
••3310  XS=''X" 

13320  YS=-5IF" 

'3330  oab$^''STRES3  INT  FACTOR  DIST." 
13348  END  IF 
13350  IF  Nd=3  then 
133G0  XS=-X" 

'  3370  yS=''RATE" 

'3330  LaD$  =  ''CRACK  GROWTH  RATES." 

13390  END  IF 

13400  CALi.  P  U  I  1  Pld'.  •  )  .Labs  ,XS  ,YS  ,Nn  ) 
13410  ' 

13420  GCLEAR 
1  3430  OUTPUT  2  i  ■  K  ; 

1 3440  GOTO  12990 
13450  SUBENC 
' j460  ' 


13470  I . . . . . 

3480  I . 

13490  I 

13500  SUB  Pit H Pt5< •  » .StrS .XlabelS .YlabelS  .Nn  ) 
13510  GINIT 
13520  graphics  on 
1  3530  xriar  =  P  t  5  (  I  ,1  .1  ) 

'  "  d  C'  y  r*  1  n  =Y  p  t  s,  i  }  I  M 
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I  35bt5  Yna>  =r  t  s  I  I  ,  i 

13560  Yinin  =  Pt5x  !  ,I  .3) 

>3570  for  Gg=i  TO  Nn 

I3SS0  FOR  Hh'0  TO  P t s i Gg  , 0  .0  )+ 1 

13590  IF  Pt5(  6g  ,Hh  .1  ) 'Xna»  THEN  X >ia>  =P t  s c  Gg  ,Hh  ,  !  > 

13600  IF  PtatGg  ,Hh  ,  1  )  .Xciin  THEN  Xr-. i n=P t  b t Gg  .Hr- , i 

I3EI0  IF  Ptst  Gg  ,Hh  ,0  ) 'Yna/i  THEN  y  na^  =P  t  s  i  Gg  .Hn  .  3  ) 

13B20  IF  Ptsi  Gg  ,Hn  ,;  >  .Yfij  n  THEN  Ynin^Pts'Gg  Hr.  O' 

13630  NEXT  Hh 

13640  next  6g 

13650  PRINT 

13660  print 

13B70  print  “X. range-  -iX-nin;-  TO  -;Xr,a- 
13680  print  "Y.RANG£=  ":Ytriini'  TO  " -.ynn- 
13690  INPUT  "X-AXIS  LABEL  SPAC  ING-* "  ,  x  spS 

13700  IF  LENCXsDSt=0  THEN  GOTO  15690 

13710  ON  ERROR  GOTO  13690 

13730  Xsoc=UALi Xsp$  1 

1 3730  OFF  ERROR 

13740  input  'Y-axis  label  spacin67- ..sns 

13750  IF  LEN<YsdS)“0  THEN  GOTO  13740 
13760  ON  ERROR  GOTO  13740 
13770  Yapc-L'ALv  Ysp*  ) 

13780  OFP  ERROR 
'3790  OUTPUT  0i"K''i 
13300  MOUE  35,95 
13810  label  Str$ 

13830  UIEUPORT  15,115,30.90 
1 p830  FRAME 
13840  X  d  1  f  f  ®Xfnax -Xfb i n 
'3650  Y  d  1  T  f  =  Y!'>ax  -  Yns  n 

13850  WINDOW  Xmin  .Xnax  +  Xpi  f  f  /  I  0  ,  Ynin-Ydi  f  f  / 1  0  Y<^ax*  Y d i  f  f  / 1  0 
13870  AXES  Xspc/S ,Yspc/5 ,0  ,0  .5  .5 , 1 
13880  CLIP  OFF 
13890  CSIZE  3 

13900  FOR  XlaB=0  TO  Xmax+X di f f / l 3  STEP  Xspc 
13910  IF  Xlabamin-Xdif  f/l3  THEN  GOTO  13940 
13930  MOUE  Xlab-Xdif f/l40,YmiP-3.Yd.ff/l0 
13930  LA3EL  USING  "S,K"ixiab 
13940  NEXT  Xlab 

'3950  MOVE  Xxiax  +  Xdi  f  f /8  ,  Yn  1  n-Yd  1  f  f /5 
13950  LABEL  USING  ‘  *  ,K  "  ;  X  labe  1* 

13970  FOP  Ylao=0  TO  Ynax  +  Ydi f f / I  3  STEP  Yspe 
13980  IF  Ylab‘,YMir-Ydiff/t3  THEN  GOTO  14010 
13930  MOUE  Xnin-Xdi f f / I  0  ,  Yiab-Ydi f f /30 

14000  label  using  “S.K'iYiab 
14010  next  Yiab 

14030  FOR  Ylab=0  TO  Yn i n-Yd j f f / 1 0  STEP  -Yspe 
14030  MOUE  Xnin-Xdif f/10  ,Ylab-Ydiff /3e 
14040  label  using  •»,K"iYiah 
14050  NEXT  Y4ab 

I40B0  MOUE  Xiiin-Xdi  f  f  /  10  ,Yr>a>  +  Ydi  f  f /8  4I 

14070  label  using  '»  ,K- sYiabel* 

14080  P-1 

14090  FOR  Oq*l  TO  Nn 
'4100  line  TYPE  Oq+3 
14110  FOR  Rr»0  TO  Pts(OQ.0,0) 

14130  PLOT  Ptst  Qg  .Rr  ,  I  )  ,Pts(  Og  .Rr  ,3  )  ,P 

14130  NEXT  Rr 

14140  MOUE  0,0 

14150  NEXT  Oq 

14  151  QS-" 

14160  INPUT  •C0PY7"_0$ 

M-’-'x  IF  OS='Y"  THEN  GOTO  142B0 
I  IF  0*='N'  THEN  GOTO  14380 

'  GOTO  14151 
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1 4250  ' 

1 4260  OUTPUT  2 i "M" ; 

14270  I 
14280  5UBEN0 
14290  ' 

14  300  . . . . . . 

14310  I . . . . . * . . . 

14320  I 

14330  SUB  Dp io t ( Gaus3< •  )  ,U1 1  < 

14331  IDELTA  K  AND  R 

14332  I 

14338  Xi>ia>=Gaus5(  0  ,  I  ) 

14339  Xr!in=Gauss<  0 , 1  > 

14540  Y  1  ma;<-Gaus5 <  0 , 4  ) 

14341  Y2min  =  GauS5< 0  ,5  ) 

14344  FOR  Pp=0  TC  Ult 

14345  IF  Gauss!  Pp  ,  I  )>Xma>  THEN  Xi<)a>  ^Gaussl  Pp  .  I  > 

14  34S  IF  GaussCPp.l)'. Xpin  THEN  Xmin^GausslPp,!) 

14347  IF  Gauss!  Pp  .4  ) -Y  1  isas  THEN  Y  I  na- “Gauss!  Pp  ,4  ) 

14348  IF  Gauss!  Pp  .5  Y2nin  THEN  y2nin“6auss!  Pp  ,5  ) 

1435!  NEXT  Pp 

14353  PRINT  •■XRANGE=  0  TO  ";Xna;> 

14354  INPUT  "X-AXIS  LABEL  5PAC ING’ " ,X spc 
1  43G2  '■ 

14363  OUTPUT  2i "K"  1 
I  4  355  G I N I T 

14366  graphics  on 

14367  UIEWPORT  15.115.50,90 

14369  I 

14370  UiiNDOW  0  ,Xnax  Y 1  na>,  ,  Y  1  maA 
1437!  AXES  X5pc/5,2.0,0,5,5,! 

1 4  3~2  AXES  X spc /5  , !  , 1  - Y2n 1 n  ) • Y 1 pa> / ■ - Y2n I n  >-Y 1 ma-  .Xnax  ,0  ,5 , 5  ,  i 

14373  CLIP  OFF 

14374  CSI2E  3 

14375  FOR  Xlab=0  TO  Xnax  STEP  Xspc 
I437G  MOUE  Xlab-.0l  ,-Y I  max -Y 1 na < / 1 0 

14377  LABEL  USING  ''t,K"iXiab 

14378  NEXT  xlab 

'4379  MOUE  X  isax-X  max /' I  0  ,  -  Y  1  nax  -  Y  1  ma  . '5 
1  4380  LABEL  USING  "S  .K"  i  "L.  " 

14381  IP  Yl»*iax<10  THEN  Ylinc“l 
14392  IP  Ylnax.'  =  l0  THEN  Y!inc=5 
14384  FOR  Yiab=0  TO  Ylnax  STEP  Ylinc 
I  43S5  N0U2  0- .  1  •  Xxiax  ,  Y I  ab- .  5 
14386  LABEL  USING  -t.K'jYlab 
143B7  NEXT  Ylab 

14388  MOUE  -Xna-./?  .Ylmax  +  YInaA /10 

14389  LABEL  USING  "  »  ,K  '  :  "  Kmax -Kn  i  n  . 

14390  ' 

14391  IF  Y2nin,.7  THEN 

14392  Y2inc“. I 

14393  ELSE 

14394  Y2inc=.S 
14335  END  IF 

14396  FOR  R2ab“0  TO  -Y2nin  STEP  YZinc 

14397  MODE  Xnax+Xnax/l5,!Riab  +  !y2nin/20)-y2riin)*Ylnax*2/(-Y2r,in»2>-Ylr)aA 

14398  LABEL  USING  -t.K-iRlab 

14399  NEXT  Rlab 

14400  FOR  Rlab“0  TO  Y2nin  STEP  -Y2inc 

14401  MOUE  Xnax  +  Xf>iax/15,!Rlab+!  Y2nin/20  )-y2riin)»Ylnax»2/(  -  Y2n  i  n»  2  )- Y  I  max 

14402  LABEL  USING  "t.K'jRlab 
I  4403  NEXT  Rlab 

14404  MOUE  Xnax +X{*iax  /  1  0  ,  Y 1  max xY  1  pidx  '  1  0 

14405  label  USING  ,K'"  :  "R.  " 

14406  NOUE  Gau5s!0,l),GauS5!0,4) 

14407  FOR  Ff=0  TO  Ult 
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I440BHLUI  uauSS'rr.ii.adUsavri.H  .  . 

14409  NEXT  Ff 

14410  MOUE  Gau55(0,l  1  ,  (  Gauss(  0  ,5  j-Y2nin  Imax/i )-Yl  i^a,- 

14411  LINE  TYPE  4 

14412  FOR  Gg=0  TO  UU 

14413  PLOT  GausstGg.l  ),(  Gauss  (  Gg  ,5  )-Y2"ii  n  )•  2  •  Y  I  na;-. /( -2  •Y2nin  ;-Y  I  na  •. 

14414  NEXT  Gg 

14415  I 

I  44  IB  Q$="" 

14417  INPUT  ■COPY'’", OS 

14418  IF  0$="Y“  THEN 

14419  OUTPUT  2; "N" j 

14420  GOTO  14425 

14421  END  IF 

14422  IF  Q$=  N"  THEN  GOTO  14425 

14423  GOTO  1 441B 
1 4425  SUBENO 
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Fig  1 


Fig  1  Typical  residual  stress  distributions  showing  partitions 
for  interpolation  routine 


TM  MS  1207 


FIfl  3 


4 


Section  a 


Opt1  =  0 
Radius 


Section  b 


♦ 


Multi-plot - 1 _ ^ 

N  To(#) 

T 

Opti 

0  -  No  residual  stress 

1  -  RSF  input  via  keyboard 

2  -  RSF  input  from  file 

? 

Max  applied  stress 
Min  applied  stress 


Opt1  =  1  Opt1  =  2(#) 

I 

Radius  Name  of  file? 

Originof  coords?  (0.1)  , 

■5)  -  Centre  of  hole 

1  -  Edge  of  hole 

Input  coords  (blank  to  end) 

Input  X 
Input  stress 

—  Change  any?  - ^  j - 

^  Plot  points 

(X-label  spacing) 

(Y-label  spacing) 


N 


No.  of  regions? 

Positions  of  boundaries? 

Re-plot  with  interpolation 
OK? 
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Fig  3  Flow  diagram  of  computer  program 
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Plot  groups  of  distributions  (Y/N) 


Name  of  file? 


Another  file?  (Y/N) 


Section  f 


Which  do  you  want  to  plot? 


0  Exit 

1  Residual  stress  fields 

2  Stress  Intensity  factors 

3  Rate  values 

I 

Plot  choice  (1,2  or  3) 

i 

0  Exit 


Fig  3  (cont)  Flow  diagram  of  computer  program 
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Section  c 


Section  d 


Corr=  0 
No  correction 


Section  e 


! 

1 _ 


fY 

Opt  2 

(•)  Crack  lengths  at  which 
SIFs  are  to  be  calculated 

I 

1  Range  of  crack  lengths 

2  Specific  crack  lengths 


Opt2  =  1 

I 

T 

Input  range 
Start  of  range 
End  of  range 
Incremental  length 


Opt2  =  2 

1 

f 

Enter  crack  length 
(blank  to  end) 


1  f 


Which  correction  routine? 


Corr 

0  No  correction  factor 

1  Newman  and  Raju  -  G2 

2  N  and  Fl/Shiv  and  Forman 


N 


Corr=1 
N&R  -  G2 


N 


Corr  =  2 
N&R/S&F 
■-  Thickness 
Aspect  ratio  (A/C) 

—  Width  correction?  (Y/N) 
Width 


f 


— ►  Printout  K  values  etc 

t 

Copy?  (Y/N) 

PlotK^axVSL(Y/N) 

Plot  AK  &  R  vs  L  (Y/N) 

t 

Calculate  rate  values?  (Y/N) 

i 

K  (c  value? 

Copy  (Y/N) 

T 

-Store  distributions  (Y/N) 


N 


Create  new  file  (Y/N) 
File  name  (Y/N) 

I 

-  Another  run.  same  dist 
(Y/N) 


To(*) 


Fig  3  (concluded) 


0 


Dimensions  in  mm 


Fig  4  Fatigue  test  specimen 
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